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1. Quantum Confinement and Resultant Quantization of Energy Eigenvalues (1) 

2. Sawtooth Superlattice with and without Strong External Reverse Bias 

a. Zero-field sawtooth superlattice 

b. Strong external reverse bias F ap , resulting in staircase superlattice 


4. Transfer Matrix at a Potential Step 

5. Transfer Matrix for a Single Step Barrier: E > V<j 

a. Transmission coefficient for uniform effective mass 

b. Transmission coefficient for different effective masses inside and outside the barrier 

6. Tunneling In a Single Step Barrier: E < Vo 

a. Transmission coefficient for uniform effective masses 

b. Transmission coefficient for different effective masses inside and outside the barrier 

7. Transfer Matrix for a Single Step Barrier: E = Vo 

8. Multiple Step Barrier (MSB): N periods 

9. Multiple Step Barrier with Three Periods 

The single resonance at Eq is split into two levels at Eq ± AE by coupling between the two wells. 
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10. The Difference Between Multiple Step Barriers and Multiple Quantum Wells 
(MQW). 
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27. Esakl-Tsu Multiple Step Barrier Geometry (11) 
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b. Applied field strength F - F ap over the length I of the MSB 
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Barrier is 10 nm wide, 0.33 eV high. Effective mass is uniform and equal to free electronic mass 
mo. 
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Barrier is 10 nm wide and 0.33eV high. Effective mass is mj n = 0.0871 inside the barrier, m^t = 
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30. T-E Data for Single Step Barriers lOnm Wide and 0.33 eV High. 

a. Superposed T-E curves for GaAs-AI(0.4)Ga(0.6)As barrier, and barrie r wit h uniform e ffective 
mass m e ff = mo everywhere. Curve 1 : GaAs-AI(0.4)Ga(0.6)As, mj n = 5.0871, rriout = 0.0636; 
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34. Total Transfer Matrix Element for Single GaAs-AI(0.5)Ga(0.5)As 
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36. Bound States in Single Quantum Well. 
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agreement with equation 139. 
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field strength; the maxima in the data correspond to quasi-bound states of the quantum well. 
Maxima shift to energies deeper in the well as the field strength increases. 

Curve 0: F ap -1.0 xl O' 2 eV/nm; Curve 1 : F ap = 1 .5 xl O' 2 eV/nm; 

Curve 2: F ap = 2.0 xIO* 2 eV/nm; Curve 3: F ap = 2.5 xIO' 2 eV/nm; 

Curve 4: F ap = 3.0 xIO’ 2 eV/nm 

40. Stark Shift of Bound State Energy Level under Uniform, Un res tr icte d Field. 

From Figure 39, the energy levels of the maxima in (F/E) are plotted as function of the applied field 
strength. The bound state energies are quadratically dependent on the applied field strength. 

41. Stark Shift of Bound State Energy Levels under Uniform, Unrestricted Field. 

From Figure 40, the Stark shift relative to the zero-field level is plotted as a function of the applied 
field strength. Also shown are the data of Austin and Jaros (41) for an identical quantum well. 

42. Stark Shifts under Localized and Unrestricted Uniform Fields Compared. 
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43. T-E Curves for Two-, Three-, and Five-Barrier GaAs-AI(0.5)Ga(0.5)As 
Multiple Step Barriers. 
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45. T-E Curves for Two-, Three-, and Five-Barrier MSB's Calculated By Esaki 
and Tsu (11). 
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46. J-V Curves for Two- and Three-Barrier GaAs-AI(0.5)Ga(0.5)As MSB's 

These curves are calculated for the same MSB's whose T-E curves appear in Figure 43. 
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47. J-V Curves for Two- and Three-Barrier MSB's Calculated by Esaki and Tsu 

( 11 ). 

Note the agreement between these curves and those of Figure 46. 
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1. INTRODUCTION 


Recently there has been much interest in semiconductor heterojunction superlattices 
both theoretically and practically. The theoretical interest fies in the fact that quantum effects are 
observable in these macroscopic structures (1-3). The practical interest lies in the application of 
these superlattices as novel electronic devices with desirable characteristics (3-8). 

The quantum mechanical effects displayed by semiconductor superlattices are 
manifestations of both quantum size effects and tunneling (1-11). Quantum confinement of 
electrons and holes in structures which are of the same order as the de Broglie wavelength of the 
particle results in quantization of the energy eigenvalues where previously a continuum existed. 
The degree of quantization depends on the number of dimensions in which the confinement 
exists. The motion of the confined particle then has a reduced dimensionality depending on 
whether the quantum confinement is present in more than one dimension. For instance, 
confinement in one direction results in approximately two-dimensional motion with partial 
quantization of energy eigenvalues. If the particle is confined in all three dimensions, the energy 
levels can be sharply quantized. (See Figure 1) 

Tunneling is the process whereby a quantum particle can cross a potential barrier, which 
classically would be completely impenetrable because of its higher energy, to a state of equal or 
lower energy. The wave function of the particle extends into or even through the barrier when the 
barrier potential is finite. A non-zero particle current density through the barrier then results from 
the tunneling process. Quantum mechanical tunneling has been the basis of many 
semiconductor devices, starting with the Esaki tunnel diode, named for its inventor Leo Esaki, 
and introduced in 1958 (12). A discussion of the evolution of tunneling theory from 1928 up to 
the early seventies has been presented by L. Esaki (3, pp. 47-77). 

When a particle interacts with and is confined by two or more barriers of finite height and 
not too great thickness, its wave is reflected multiply off each potential barrier reached by 
tunneling. When the confining region’s dimension is some multiple of the wavelength, the particle 
"resonates" in the regions where its energy is greater than the local potential. At these 
wavelengths the tunneling current is amplified. Actually, both size quantization effects and 
resonance result from the same source: the constructive interference of forward and backward 
waves. This is the source of the quantized energy levels that result from the confinement of the 
particle. 

Resonant tunneling figures prominently in the transport of carriers through 
semiconductor superlattices, and to understand it is not only desirable theoretically, but is also 
central to the application of these structures as electronic devices. Resonant tunneling of 
electrons and holes in the conduction and valence bands leads to formation of sub- or mini- 
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bands, whose widths and energy levels will affect the operation of superlattice devices. Tunneling 
calculations for semiconductor superlattices can provide this information, which can then be 
implemented in design of improved structures. By finding the transmission coefficient through 
the structure as a function of energy for either carrier, one can locate the energy levels of the 
resonances. It is also possible to find the effects of an external electric field on the transmission 
coefficient, and hence upon the resonance energy levels. Much effort has already been directed 
to this end for step superlattices; however, the sawtooth or graded band-gap superlattice 
proposed by F. Capasso (3) has up until now not been studied theoretically, except for an analysis 
of the multiplication noise associated with its use as a photo-detector. This is the aim of the 
present study. 

A resurgence of interest in resonant tunneling in heterojunction semiconductors has 
been spurred by recent advances in molecular beam epitaxy, ( MBE), w hich provides abrupt 
interfaces (on the order of a monolayer) as well as very uniform layer thicknesses (4-8,13). MBE is 
used to make superlattices of multiple barriers and wells in which the energy levels of the virtual 
states are consistent from well to well, leading to miniband formation and therefore to efficient 
transport of carriers through the structure. The energy levels in quantum wells and superlattices 
have been most recently and thoroughly reviewed by Altarelli (2). Excellent discussions of the 
devices which can be made from semiconductor superlattices grown by MBE, and the relationship 
of the superlattice structure to energy levels and resultant device performance are given by 
Capasso et alia (4-8). 

The earliest theoretical exploration of resonant tunneling in semiconductor 
heterojunction superlatiices was published by L. Esaki and R. Tsu in 1970 (9), who predicted 
negative conductance, caused by electron transport into negative effective mass regions of the 
minizone, and Bloch oscillations. The first experimental observation of resonant tunneling 
through a double barrier was made by Chang, Esaki, and Tsu in 1974 (10). 

A fairly comprehensive body of work addressing carrier transport through GaAs-AIGaAs 
step superlattices exists with which comparisons may be made when studying other kinds of 
superlattice, starting with the work of Esaki and Tsu cited above. In a paper published in 1973 
(It), they presented the calculated transmission coefficient as a function of incident energy, and 
tunneling current as a function of applied voltage, for two, three, and five GaAs-AI(0.5)Ga(0.5)As 
step barriers. A transfer matrix method was used for this work. Their computations were for thin 
layers (50 and 20 A wells and barriers, respectively), neglecting the potential gradients caused by 
the applied field within each layer, but taking account of the total potential drop between 
successive periods. They took account of effective mass variations throughout the superlattice. 
The correlation of the current peaks in the calculated J-V curves with the resonance energy levels 
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in the transmission data indicates a Stark shift in the levels that is linear in the electric field 
strength. 

More recently, Marsh (14) applied an empirical pseudopotential formulation to tunneling 
calculations through a GaAs-AI(x)Ga(1-x)As double heterostructure step barrier under zero 
external electric field, and compared the results to those obtained by the effective mass 
approximation. In taking account of the effective mass variations throughout the heterostructure, 
Marsh made the distinction between an effective mass approximation which uses conventional 
wave function connection rules at interfaces where there is a sharp change in effective mass, and 
one which uses the modified connection rules suggested by Kroemer and Zhu (15,16). The 
latter, about which more will be said below, maintain particle flux continuity through the interface. 
Marsh found good agreement between the effective mass approximation and the empirical 
pseudopotential when the modified connection rules are used, and the aluminum concentration 
in the barrier layer is such that it is still a direct gap material. Significant discrepancies were 
observed between the two methods when direct-indirect interfaces were studied, and in all cases 
when the unmodified connection rules were used. 

A modified scattering matrix formalism was recently applied to calculating the resonance 
energy spectra for electrons in multiple GaAs-AIGaAs quantum wells (17). The Stark shift of the 
levels was also calculated and found to be linear in the field strength. Although in that work the 
effective masses appropriate to each region were used, it appears that the conventional wave- 
function connection rules were used in deriving the scattering matrices, which may have resulted 
in some error in the computed energy levels. 

MBE growth of these structures provides very precise control of the aluminum 
concentration in AIGaAs even at the monolayer level, allowing linear and even parabolic grading of 
the band-gap. Several novel electronic devices using both linearly and parabolically graded-gap 
materials have recently been described in the literature (4-6). A superlattice proposed by Capasso 
for use as a low-noise, high-gain solid-state avalanche photodiode is described in (7). In this 
structure, the conduction and valence band edges in each stage describe a sawtooth profile, that 
is, the band-gap in each layer is approximately linearly graded, (refer to Figure 2) This is achieved 
by varying the aluminum content linearly within the layer.At each interface there is a band-gap 
mismatch between pure GaAs and AI(x)Ga(1 -x) As which is taken up mostly in the alignment of the 
conduction bands. Five or six stages are contained between n and p doped GaAs cladding layers. 
In operation, photo-electrons are accelerated down the structure by a strong external reverse 
bias. At each interface, because of the built-in potential drop afforded by the conduction band 
misalignment, the electron acquires enough kinetic energy to impact-ionize a lattice atom and 
liberate another electron. In this way the electrons multiply at each stage. The small valence band 
discontinuity in AIGaAs prevents the impact ionization of the holes. The larger the conduction 
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band discontinuity, the higher the probability that each electron will impact-ionize at each stage. 
This means that very low-noise multiplication can be achieved, approaching the light- detection 
performance of a photomultiplier tube. 

The size of the potential mismatch at each interface in this device is determined by the 
amount of aluminum in the AI(xjGa(1-x)As at the end of each stage. For x greater than 
approximately 0.4 (18,19), AI(x)Ga(1-x)As is an indirect gap material; pure AlAs has an indirect 
bandgap, and pure GaAs is a direct-gap material. The problem of electron tunneling through a 
direct- indirect gap interface must be handled differently than the simpler case of transport 
between two direct-gap layers, as the work by Marsh cited above indicates. Physically, the 
transport of electrons may be hindered by the competition of the direct and indirect bands in the 
high aluminum material. This may partially undo the advantage of the large conduction band 
discontinuity associated with high aluminum content. Hence consideration of AI(x)Ga(1-x)As 
sawtooth superlattices was limited in this study to compositions with x no greater than 0.4. 

Sawtooth superlattices achieved by compositional grading differ from step superlattices in 
several important ways. One effect of the compositional grading is that the effective mass of the 
carriers is nowhere constant throughout the structure. In step superlattices the effective mass is 
constant within each layer, but there is a sharp change at each boundary. In sawtooth 
superlattices in AIGaAs the effective mass of the electron in the direct band is a quadratic function 
of the amount of aluminum. Since the aluminum content is linearly graded, there is a quadratic 
variation of the effective mass and the band-gap within each layer, and a large discontinuity at the 
end of each stage as well (18,19). 

Another significant difference between the step and sawtooth superlattices is the 
presence of internal or "quasi" fields associated with the band-gap grading. One of these is a 
constant field which is the gradient of the band-gap. This field accelerates electrons and holes in 
the direction of the narrowest bandgap, and opposes the reverse bias applied to the device. 


f = 



(la) 


F 


h 



(1b) 


Ec(z) and Ev(z) are the conduction and valence band edge energies, respectively. This field is 
practically much less than the applied reverse bias for direct-gap AIGaAs-GaAs sawtooth 
superlattices. 
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Furthermore, the varying effective mass of the carriers gives rise to another quasi-field 
which assists the external reverse bias in accelerating carriers through the structure. This field is 
given by: 


F e = 


F h = 



(2a) 


(2b) 


where m eeW (z) and mhen(z) are the electron and hole effective masses, respectively. For electrons 
in a direct-gap AIGaAs-GaAs sawtooth, this field strength is much less than that due to the band- 
gap gradient, and these two quasi-fields offset one another somewhat. 

One more major difference between step and sawtooth superlattices is the way they 
deform under high external fields. Figure 3 illustrates this difference. In a step superlattice, no 
matter how strong the field, the barriers remain such that there is always a well between them. 
Resonances can in principal always be formed in these wells even under strong fields. The 
sawtooth, on the other hand, becomes a staircase structure at applied voltages greater than the 
sum of the conduction band discontinuities over all the stages. Then the triangular wells between 
the sawtooth barriers no longer exist. The electron still interacts with the staircase, but tunneling is 
no longer occurring between the barriers. Therefore it is expected that the transition from 
sawtooth to staircase should be signalled by some feature in the cunrent-voltage characteristic. 

An analysis of tunneling through a graded gap superlattice should thus take account of 
not only the external electric field and structural parameters such as number of stages, layer 
widths, composition, and interfacial conduction band discontinuities, but also of the quasi-fields 
caused by the gradients in the band-gap and effective mass. The differences between step and 
sawtooth superlattices should be reflected in the results of tunneling calculations. 

It is the central purpose of the work described here to study the resonant energy levels of 
graded band gap GaAs-AI(x)Ga(1-x)As superlattices for x less than or equal to 0.4 (i.e. for direct 
gap material only), as a function of the applied and quasi-fields discussed above. A transfer matrix 
method is used to calculate the transmission coefficient and tunneling current for conduction 
electrons in sawtooth superlattices similar to those which might be used as avalanche 
photodiodes. Tunneling calculations are also presented for the same kind of step superlattice 
investigated by Tsu and Esaki (11). These are shown to agree with those previously obtained, 
and are also used to illustrate the differences and similarities between sawtooth and step 
superlattices with similar compositions and structural parameters. Throughout this work particular 
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attention is paid to role of the effective mass in electron transport through these structures. The 
justification for using the modified wave function connection rules when effective mass 
discontinuities are encountered at interfaces, is discussed. The effects of neglecting effective 
mass variations, or using inappropriate connection rules, are also explored. 
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2. THEORY OF TUNNELING IN SUPERLATTICES 


2. A. Quantum Mechanical Tunneling 

In quantum mechanical tunneling, a particle of energy E is incident on one or more 
potential barriers of arbitrary shape and height.lts behavior is described by the Schroedinger 
equation; specifically, the wave function inside and outside the barrier, energy eigenvalues and 
resonances, the transmission and reflection coefficients, and transmitted and reflected probability 
density currents may be calculated. In the problem of tunneling in superlattices, the salient 
features are essentially one dimensional: most often the superlattice is formed of parallel layers 
alternating in one dimension, say the z direction. The tunneling barriers then extend in the z 
direction, while momenta in the x and y directions are constants of the motion. The one- 
dimensional Schroedinger equation is then used to describe the motion of the tunneling particle: 

#L± + (V(z) - E) 4> = o (3) 

2m eff dz 2 

where y is the wave function, V(z) is the potential due to the superlattice, and m e ff is the mass of 
the incident particle (20,21). 

In superlattices formed of real solids, the mass in the expression above is the effective 
mass appropriate to each layer. In the effective mass approximation, the ionic potential of the 
crystalline lattice is not dealt with directly, but instead is taken account of by the parameters of the 
effective mass and the energy band edge.(22) 

The solution to the Schroedinger equation in regions of constant potential where V(z) = 
Vj, is the set of plane waves: 

ik:Z — ik ; z 

i//j= Aje 1 + Bjfl * (4) 


where; 


_ V^.w'E-Vj) 15) 

j t. 

When the energy E is greater than the potential Vj, (E - Vj) > 0 and the general solution 
above is composed of plane waves propagating to the right: exp(ikj z), and to the left: exp(-ikj z). 
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The particle propagates throughout this region like a free particle. This is characteristic of the 
motion of the particle in the regions outside the barriers. 

When the particle encounters regions where its energy is lower than the potential, then 
(E - Vj) < 0, and the solution to the Schroedinger equation is the sum of exponentially growing 

and decaying parts: 


where Kj is: 


k i = ik i z 


i N /2m eff (E-V j ) 

ft 


(7) 


This wave function describes the penetration of the particle into a barrier. When the effective 
mass m e ff of the particle is constant from region to region, the coefficients Aj and Bj are found by 
applying the standard connection rules of matching the wave functions and first derivatives at 
each boundary: 


l^(z) = 

*j+1 (z) 

(8a) 

*•(*> = 


(8b) 


In semiconductors the effective mass is dependent on composition. In semiconductor 
superlattices, the effective mass therefore makes a discontinuous jump at each interface between 
dissimilar materials. When this is the case, the standard wave function matching procedure above 
will not result in conservation of the probability density current through the interface (15,16,23). 
The connection rules must be modified so that this current is conserved. One approach (15,16) to 
this is to redefine the wave functions on either side of the interface as: 



Once this is done, the standard practice of matching the wave functions and their first derivatives 
can be applied to the Xj to obtain the coefficients Aj and Bj. To summarize, the former wave 
function y is not continuous at each boundary when the effective mass of the particle changes 
abruptly there, but the renormalized wave function x. used with the standard connection rules, 
results in coefficients which maintain the continuity of particle flux J through the boundaries: 


8 


i will 1 1 1 


J i - Vi 


( 10 ) 


-it) / * dXj dx* \ -ij> / * dx j+1 d X j+l \ 

dz ' * dT) ’ 2m 0 \Vl ~ Vi 1‘ I 

2.B. Transfer matrix method 


(ID 


2.B.I. General properties of the transfer matrix 


Transfer matrices can be derived describing the propagation of a wave through a 
superlattice of wells and barriers. The formalism is identical to that used in optics, where ray 
matrices are applied to the propagation of light through optical elements. (24) 

The transfer matrix M used in the present work is closely related to the scattering matrix S 
that is applied to three-dimensional problems of nuclear scattering. The S matrix is most useful for 
formulating symmetry properties, whereas the M matrix Is best applied to the one-dimensional 
problems dealt with here. It has been used in tunneling calculations for semiconductor multiple 
step barriers (MSB's), by Esaki and Tsu (11) and C. Schwartz (17), and for electrode-polymer 
interfacial layers by Meijer and Van Roggen (25,26). 

The wave functions on either side of a potential step at z = 0 are written: 


V 

i k:Z — 

Aj e 1 + Bj e 

i kjZ 

z < 0 

(12a) 

Vl 

ik i+1 z 

= A j+1 e i +1 + 

B e _ik i +lZ 
B j+1 e 

N 

V 

o 

(12b) 


There are two linear, homogeneous equations relating the coefficients on either side of the step. 
The transfer matrix M expresses these equations: 



(13) 


The origin of the My are the matching conditions at the step. M is determined to be: 



( 14 ) 
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Note that mj/mj + i is a factor of unity when there is no effective mass change at the 
interface. An overall transfer matrix through an arbitrary number N of steps, is just the product of 
the N such matrices found by matching at each step: 

N 

M = n M: (15) 

j = i 1 


Some general statements may be made about the properties of M, without knowing the 
details of its form, as long as the potential to which it is applied is real, goes to zero at infinity, and is 
symmetric about the origin. {See reference 20 for the discussion to follow.) It i s importa nt to note 
that the latter two conditions do not apply when a uniform electric field is present, nor when there 
is an asymmetric effective mass gradient throughout the structure, as is the case for t he sawtooth 
superlattice. In that case, and in the case of the MSB under a uniform electric field, tunneling 
calculations are complicated somewhat by this lack of symmetry. 

The scattering matrix S for a one-dimensional potential is: 



This matrix relates the magnitudes of the outgoing waves to those of the incoming ones, and is 
easily related to the M matrix when the probability current density is conserved. In that case, S can 
be shown to be unitary, and symmetric as a consequence of the time reversibility of the 
Schroedinger equation for a real potential. M is then related to S as: 



When the potential is even, the Schroedinger equation is invariant under space 
reflection. Then Si i = S22 and Si 2 = S21 . This fact combined with the unitarity and symmetry of S, 
leads to the relationships between the My: 
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M 11 - M *22 


m 12 = M* 12 = -M 12 = M* 
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det (M) = 1 

To summarize, the expressions above apply to transfer matrices for real, even, potentials which go 
to zero at infinity, and for which the probability current density is conserved. 

Making the assumption that there is no wave incident on the barrier from the right implies 
that Bj + i = 0. Then the transmission coefficient T, defined as the magnitude squared of the ratio 
of the transmitted to the incident wave amplitudes Aj + i/Aj, is: 


t = 





( 18 ) 


When Mu * 1. then the barrier is effectively transparent to the incident particle, and 
transmission is at a maximum of 1. This phenomenon is known as resonance. When Mu is very 
large, the barrier is effectively opaque, and most of the incident energy is reflected. Thus the 
transmission coefficient is related to the total transfer matrix M via the element Mi i . 

The transfer matrix is also used to calculate the ratio of the transmitted to the incident 
current densities. Applying equation 11 for the probability current density, to the incident and 
transmitted components of the wave functions Ai and At on either side of the interface, gives the 
incident and transmitted current densities: 

j= = | a. | 2 (19a) 

m i 

J, • ^|A ,| 2 »»> 

m t 

The ratio of the transmitted to incident current density is then: 

J is often referred to as the tunneling current. It is the actual property of physical interest in 
tunneling problems, and is proportional to the transmission coefficient. Thus the element Mi i of 
the total transfer matrix determines the tunneling current through semiconductor superlattices. 

The presence of the effective mass should be noted in the expression above for the 
tunneling current: the ratio of the effective masses in the initial and final media will affect the 
current. Tunneling from high to low effective mass is less probable than tunneling in the opposite 
direction, for instance. This is anticipated in equation 2 for the quasi-field caused by effective 
mass variations: 
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v .i kT a, 

2 \ % 

The potential difference between the initial and final media due to the effective mass gradient is 
given by: 

AV = j 

Inspection of the equation above shows that electrons will be accelerated toward the medium with 
the higher effective mass. 

2.B.2. Tunneling: Single Step Barrier 



The single step barrier is the building block for the most common kind of superlattice, the 
multiple step barrier (MSB), often referred to in the literature as the multiple quantum well (MQW). 
In the absence of an electric field, the potential is constant inside and outside the barrier, and 
simple analytical expressions are obtainable for the transmission coefficient. Even when there is 
an effective mass difference between the barrier and its surroundings, this only adds a small 
constant term to the potential in each region, and analytical expressions for T are still easily found. 
To illustrate the use of the transfer matrix in tunneling problems, the transmission coefficient is 
found below for the step barrier. This will provide a basis for comparison of tunneling through 
single and multiple barriers. (See Figure 5) 

The total transfer matrix through the barrier is the product of the transfer matrices at each 

step: 

M = M 1 • M 2 (22> 


where Mi is evaluated at z = zi and is: 



( 23 ) 


and: 
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M2 is evaluated at z = Z2: 



k 2 = 


>/2m 2 (E - V Q ) 


t. 


M 


2 



— k 2^*2 


+ k 2 )z 2 



( 24 ) 


2.B.2.a. Single Step Barrier: E > Vq 


Consider first the case where the incident energy E is greater than the barrier height, E > 
Vo- (Refer to Figure 5 a) In this case ki and k2 are both real and when Mi and M2 are multiplied, the 
matrix element Mu is: 


M 
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/ 2 ( k 1 + k 2 ) 

\ 2 k 1 k 2 


sin 



(25) 


where a is equal to the barrier width (Z2-Z1), and |Mi 1 1 is, after a little rearrangement: 


M 
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4 k i k 2 + (ki| — k^j^sin^kga 


(26) 


4k^ k| 


The transmission coefficient T is thus 1/|M-j i| or: 


T = 1 + 


k^ — k2j sin* k2a 
4k^ k| 


r 


(27) 


Inspection shows that T is at its maximum of 1 whenever sin(k2a) is zero, or when k2a = rot. This is 
resonance, in which standing waves exist in the region of the barrier. At resonance, an incident 
wave packet spends a comparatively long time in the vicinity of the barrier, even though its energy 
is greater than the barrier height. 

T can be expressed explicitly as a function of incident energy E if the definitions of the 
wavevectors ki and k2 are used in the expression above. At this point the effective masses 
should be retained in the expressions for the two wavevectors: if the effective masses are 
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different inside and outside the barrier, T(E) will differ from the standard expression which 
assumes that the mass is everywhere constant. The expression for T(E) which takes account of 
the effective mass variations is: 


T = 


i + 


( E(m^ - rr^) + rr^V Q j 2 sin 2 


-1 


(28) 


4m.jm2 E(E - V 0 ) 


Note that when the effective mass is equal throughout, this expression for T(E) becomes 
the standard formula: 


T 


1 + 


V 2 sin 2 


(fs/ 


-1 


4E (E 


2m 2 (E -V Q ) 


(29) 


The difference between these two equations is worth noting. A term involving the 
incident energy and the difference between the two effective masses survives in the 
denominator of the first equation. This term changes the minimum values of the transmission 
coefficient relative to those in the constant mass case: the relative magnitudes of mi and m 2 
determine whether the minima are raised or lowered. It is also evident that the oscillatory term in 
the denominator is a function of the effective mass inside the barrier, m 2 . 

In the constant effective mass expression, the energies at which the transmission 
coefficient is maximized are easily obtained by setting the argument of the sine equal to nn: 



Inspection shows that the first maximum in T occurs for E greater than the barrier height by 
hrt/2ma . Successive maxima are scaled as n^. However, when the masses inside and outside the 
barrier are different, the location of the transmission maxima are obtained by minimizing the 
denominator in equation 29. These maxima or resonance peaks will generally not occur at the 
same energies as in the constant-mass case. Effective mass variations therefore can significantly 
affect the transmission coefficient for single step barriers even when the incident energy is well 
above the barrier height. 
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2.B.2.b. Single Step Barrier: E < Vq 


Tunneling may occur when the incident energy E is less than the barrier height Vo- (Refer 
to Figure 6b) In this case the wavevector k 2 is imaginary, and equal to: 

. _ . V / 2m eff (V 0 -E) 

i<2 i ^ * i' 

Substituting iK for k 2 in equation 28, and using the definitions of k^ and k results in the 
expression for the tunneling transmission coefficient through the barrier: 

/ ( E(m 1 - m 2 ) + m 2 V 0 sin h 2 yfi m 2 (V Q - E)) 

T = 1 + J L 

\ 4m 1 m 2 E(V q — E) 



Just as in the preceding case where E > Vq, the effective masses inside and outside the 
barrier have been retained. If the effective mass is constant, then the familiar expression for T(E) is 
retrieved: 


T 



sinh' 


2 (U 


2m «ff ,V o - 


4E (V Q - E) 



(32) 


Again, both the absolute values and the relative differences between the effective 
masses affect the transmission coefficient. In both expressions, T increases monotonically from 
zero, for incident energy of zero, to its maximum value (for energies less than the barrier height) 
when E * Vq. However, the rate at which T increases is a function of m-| and m 2 . The maximum 
value can be found by taking the limit in equation 32 as (E - Vo) approaches zero: 

T - (,<• 1331 
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It is interesting that the maximum value of T in the limit as E approaches Vq, is a function of the 
effective mass outside the barrier, and not the mass inside the barrier. 

2.B.2.C. Single Step Barrier: Lim E->Vo 

It is worth examining what happens when the particle energy coincides exactly with the 
barrier height, a condition shown in Figure 7, and demonstrating that this does not result in any 
discontinuity in the transmission coefficient, since it is a frequent occurrence when an arbitrary 
potential is broken down into many small potential steps. This is the classical "turning point" of a 
particle in a potential. It is most easily treated by going back to the Schroedinger equation, and 
immediately using the fact that E = Vo- What results is the Laplace equation over the region of the 
barrier: 

2 

=0 0 < z < a (34) 

dz 2 

The solutions to the Schroedinger equation outside the barrier are still the plane wave solutions 
described above. The wave functions inside and outside the barrier are then: 

(35a) 
(35b) 
(35c) 

(36) 

a gives the matrices 


( 37 ) 

and: 



= C + Dz 0 < z < a 

i k., z — i k, z 

\p = Fe 1 + Ge 1 z > a 

The transfer matrix through the barrier is still given by: 

C)- 

where matching the wave functions and their first derivatives at z = 0 and z = 
Mi and M 2 : 


F 
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1 / m~ 

M 2 = — / ~ L 
2 2 V m, 


1(1 - i k 1 z 2 ) e' kl 2 (l+ik 1 z 2 )e 1 ^ 


i k-j e 


ik 1 z 2 


— i k-j e 


-ik 


1 z 2 


(38) 


I * 2 = a 


From the product Mi M2, the overall transfer matrix element Mi 1 is: 

1 / i k 1 a — i k-i a \ 

=— ( {1 — i k^ a) e + e j 

The transmission coefficient is |Mi 1 |' 2 or: 


(39) 


T 



2 

m 1 E a 
2 li 2 




^iir’ 

2 fc 2 / 


(40) 


This is the same result obtained above, in the case where E < V<), by taking the limit of T as E 
approaches Vo- There is no discontinuity in T when the incident energy happens to coincide with 
barrier height, and this applies as well when the barrier consists of many small steps of constant 
potential, Vj. 


2.B.2.d. Tunneling Current Through Single Step Barrier 


The tunneling current is easily found from the transmission coefficient whether the 
incident energy is above or below the barrier height. Recalling that the tunneling current is given 
by: 



It is apparent that since the effective masses are the same on either side of the barrier, then so are 
the initial and final wavevectors. J is then just equal to the transmission coefficient. In a step barrier 
surrounded by media with different effective masses, the tunneling current would not take such a 
simple form, but would reflect the effects of the quasi-field caused by the net effective mass 
gradient. 


2.B.3. Multiple Step Barriers 
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Just as the transmission coefficient through a single step barrier is calculated, one can 
perform the same calculation for an array of step barriers or superlattice. Tunneling through an 
array of step barriers is significantly different than for a single barrier. The most striking difference is 
that resonances are possible for energies below the barrier height, whereas the transmission 
coefficient through a single barrier can only have resonances when the incident energy is at least 
equal to the barrier height. If the MSB is symmetric with respect to the origin, the transmission 
coefficient at the resonances is equal to 1 , as will be shown below, by using some of the general 
properties of the transfer matrix outlined in section 2.B.1 . 

The goal of tunneling calculations for MSB's is often to find the resonances and their 
Stark shift. It is possible to solve directly for the resonances when the MSB is not deformed by an 
external field, i.e., when there is even symmetry. When the symmetry is broken by a uniform field, 
and the Stark shift of the resonances and tunneling current density are desired, then a numerical 
solution by means of transfer matrices is useful. That problem will be dealt with in section 2.B.6; 
the solution for the zero-field case is given below for an MSB of N step barriers. 

A typical MSB is shown in Figure 8. A total transfer matrix M for the MSB can be expressed 
as the matrix Mb for a single barrier of width a, multiplied by a phase factor accounting for transfer 
over the distance w between neighboring barriers, raised to the N th power where N is the number 
of periods in the MSB: 



Sylvester's theorem (27) gives the Nth product of a square matrix as: 

n II (Mjj — Xj i) 

f(M b ) = (M b ) N = J) f(X k ) litk (43) 

k = 1 II (X k — Xj) 

i*k 


18 


■i'll ill I II Hill 111 liiUilill iiHUil | ill I III j I j. lliliiilil.il 



where Xj and X*< are distinct eigenvalues of Mb- Representing Mb as: 


■>■(: :) 


( 44 ) 


the eigenvalues are: 


A 1 2 = (u 4 - v) 


(45) 


where: 


(A + B) 


(46a) 


v J {A ~ B>: 


+ CD 


The element (Mb N )n is then, by Sylvester's theorem: 


(46b) 


iu ( A. i — B) — (X 2 — B)^ Xo^ 

' M b N| n ■ — — r— r — 

A 2 _A 1 


(47) 


and the off-diagonal element (Mb N ) 12 is: 


(M, 


(>." - »,") 


N ) = 

b 12 X -X 

x 1 -x 2 


(48) 


One can simplify these expressions by expanding X n i. 2 in binomial series and separating 
the even and odd powers into two series: 

[N+1 ] 

[N/2] / N \ 

m / V iu 


1,2 


[N/2] / N \ 2 / N \ 

(u + V) N = 2( )“ N_2q V 2q ±V 2( )u N - 2q+1 V 

q=0 \ 2q / q=1 '2q-1/ 


2q-2 (49) 
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These expansions can be used in solving for the resonances of the MSB. It is more convenient to 
use the simpler expression for the off-diagonal element, instead of (Mb N )n- 

Recall that the origin of the criterion for resonance is in the boundary condition on the 

wave: 



At resonance the coefficient b is zero because there is no reflected wave. This implies not only 
that (M t )n is 1, but also that (M t ) 12 , and by symmetry, (M t ) 2 i are zero. Then using equation 48 
above for (Mb N )i 2 gives the resonance condition for the MSB in terms of the off-diagonal element: 


M. 


12 


= 0 


<M b N ) 


12 



( 50 ) 


Hence: 


(M b N )i 2 



( 51 ) 


Using the series expansions of equation 49 for X-|N and Qives the resonance condition: 


(M b N )i2 


N+1 

2 

c s 

q=1 


| N yN-2q+1 v 2q-2 


( 52 ) 


The actual values of u and v are determined by the incident energy and the structure of the MSB. 

As an example, consider the solution for the resonances of a three period MSB such a$ 
shown in Figure 9. The expression for (Mb N )i 2 is: 


(M b 3 ) 


12 


= C 




2q —2 


( 53 ) 
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and the resonances are found by solving: 


(M h J > 


12 


= C (3 u 2 + V 2 ) 


(54) 


The fact that the determinant of Mb is equal to 1 can be used to put v in terms of u. Using a general 
expression for Mb gives: 


det (M b ) = 


A c 

D B 


(AB —) CD) = 1 


(55) 


Recalling that: 


(A+B) 

u = 

2 


V 


i 


(A - Br 


+ CD 


then: 


(A + B) 2 
4 


(A - B) 2 

+ CD = (AB-CD) = 1 

4 


or: 


u 2 — v 2 


(56) 


so that: 


(M k j > 


12 


= C(4u 2 - 1) 


u = ± 



(57) 

(58) 


In other words, there are two resonances given by u = 1/2 and u = -1/2. These are 
symmetrically spaced with respect to a single resonance in a single well between two such 
barriers. That single level is split by the coupling between wells. 

This method can be applied whenever the MSB is symmetric with respect to the origin. 
The solution for the resonances, their Stark shift, and the tunneling current when the MSB is 
distorted by a uniform field is treated in a subsequent section. 
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2.B.4. Single step potential wells: zero external field 

It is apparent that the transfer matrix method applies equally well to problems involving 
potential wells. Potential wells are taken in this study to be strictly below potential zero, or ground 
level, and a distinction is made between them and the kind of well which is formed between two 
potential step barriers.(See Figure 10) In the present work the terms "quantum well" and "multiple 
quantum well" (MQW) refer to structures below ground level. This use is adopted to avoid 
confusion, since the computation of bound states in a single finite well is different from that of the 
resonances formed between two finite step barriers. The use of the transfer matrix in finding the 
bound states of a particle in a quantum well, and the transmission coefficient of a particle 
overflying the well, is outlined below, for zero external field. 

2.B.4.a. Single Quantum Well: E > 0 (Overflying the Well) 

The treatment of a particle overflying a square potential well is very similar to the method 
already outlined for a step barrier. (Refer to Figure 11a) Given a well of width a, and depth -Vq, 
transfer matrices are found at each of its edges, and multiplied to find the overall transfer matrix M. 
The sole difference is that the wavevector in the region of the well, k 2 , becomes: 

V 2m 2< E + y (59) 

k 2 = 

The wavevector k^ outside the well is the same as for the step barrier. The transmission 
coefficient can be found by substituting k 2 above into the expression derived for the step barrier 
(equation 27). If the effective masses m^ and m 2 inside and outside the well are retained, T is 
then: 


T =11+ -5 — 


m 2 ) - m o V o) 2 sin2 




4m^ m 2 E (E + V Q ) 




(60) 


The only difference between this and equation 28 is in the algebraic sign of Vo- When the 
effective mass is constant everywhere, T becomes: 


T _ + m eff V o 2 * in2 (£V 2 m «ff< E + V o>) \ 1 

\ 4E (E + V Q ) / 


( 61 ) 
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This is also identical to equation 29 except for the sign of Vq. For 'a particle overflying a well, the 
transmission coefficient oscillates between a maximum of one and some minimum value, just as in 
the case of a particle overflying a potential barrier. Using equation 61 above with the fact that 
maxima occur when the argument of the sine is equal to me, gives the maxima in T(E): 

ti 2 n 2 n 2 I | 

E = - I V D | (62) 

2m eff a 

The physical difference between overflying a well and a barrier, is that the particle travels faster 
over the well than in the surrounding regions, and more slowly over the barrier. 

2.B.4.b. Single Quantum Well: Vo < E < 0 (Bound States) 

For a particle whose energy is between ground level and the bottom of the well, there is 
no transmission. The wave function decays exponentially into the surrounding media on both 
sides of the well, causing the energy of the particle inside the well to become quantized. The 
particle is said to be bound, and the eigenfunctions are bound states. (Refer to Figure 11b) 

To find the bound state energy eigenvalues, the total transfer matrix is set up as for the 
particle overflying the well, only now the wavevectors outside and inside the well become: 


k 


i 


- s /2m 1 (— E) 

l 


(63a) 



It is convenient to set i ki equal to k. The total transfer matrix is then: 


M well 



(k, 2 - k 2 ) 

— * sin k 0 a + 2 cos k 7 a 

k 2 k £ * 


M 
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(63b) 


(64) 


The eigenvalue condition follows from the constraint that the wave function be bounded on both 
sides of the well, and is that: 
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0 


(65) 


M 


11 


Thus the bound states inside the well are found by solving the equation: 


(^2^ k ) sin k?a + 2 cos k-. 


e ka j V*2 2. I s ' n k 2 a + 2 cos k 2 a j = 0 


( 66 ) 


k2 k 


This is a transcendental equation which can be solved semigraphically or by other numerical 
methods. 

It is now apparent why the quantum well is much different than a well formed by a couple 
of step barriers. Transmission can always occur for any particle with energy below the barrier 
height in the latter case, provided the barriers are not extremely thick, so that there are no bound 
states for such a structure. At resonance the particle is equally likely to be found anywhere inside 
the well or out of it, and off resonance, the particle is reflected from the barriers to an extent given 
by the reflection coefficient R=(1 - T). No eigenvalue condition exists such as for the quantum 
well, since any energy is actually allowed to the particle. 

2.B.5. Arbitrary potentials 


Up to this point, constant step potentials have been considered. The transfer matrix 
method can be applied to potentials of arbitrary spatial variation, however. The formalism of the 
transfer matrix requires only that the potential V(z) be a real function of z. If the potential happens 
to be symmetric then the transfer matrix has some convenient symmetries, but a lack of symmetry 
does not hinder the general approach. In fact, even for asymmetric barriers the transmission and 
reflection coefficients are the same whether the particle approaches the barrier from the right or 
from the left (21). 

In practice, when the potential V(z) is some arbitrary function of z, real but not even 
necessarily continuous, it can be approximated by steps which are piecewise constant, in the 
same manner as for numerical integration of a curve (Refer to Figure 12). Then the solutions to the 
Schroedinger equation in each interval are the simple plane waves. The functional forms of V(z) 
and the effective mass, m e ff(z), are needed so that the wavevectors to the right and left of each 

step can be calculated. Then the transfer matrix at each step is calculated from equation 14, and 
the overall transfer matrix is found for the structure by multiplication, as in equation 15. 

This method has already been used by Meijer and Van Roggen to calculate J-V curves for 
polymeric electronic devices (25,26). In such devices the electrode-polymer interfacial barriers 
and the applied electric field result in potentials which depart considerably from the simple 
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textbook step potentials discussed so far. The overall transfer matrices calculated for these 
structures yielded J-V characteristics in agreement with those which were experimentally 
observed. 

2.B.6. Field effects: Constant Electric Field 

A very important special case of an arbitrary potential is that due to a one-dimensional 
constant field. Such a potential is linear in the spatial coordinate: 

v(z) = c F ap z + v o (67) 

where F is the field strength, V 0 is a constant potential offset, and c is a constant, which equals the 
electronic charge when F is the electric field strength. This potential also describes a constant 
gravitational field, where F is the acceleration due to gravity, g, and c is the mass of a particle in the 
field. (28-30) 

Figure 13 shows a single quantum well and step barrier under a uniform electric field. 
Under such an electric field, the Schroedinger equation becomes: 


-t» 2 d 2 \p 
2 m eff dz 2 


< V o + cF ap z) * 


E* 


( 68 ) 


The solutions of this equation yield the bound states of a quantum well, and the transmission 
coefficient for a step barrier or series of barriers in the presence of a constant field. The simple 
analytical solutions obtained in the preceding sections for zero field are no longer valid. The 
transfer matrix method can still be used, but the approach adopted depends on the geometry of 
the problem: a distinction must be made between a constant field of infinite extent and one which 
is "chopped" or restricted to a certain region. In the latter case, the transfer matrix method can be 
applied in a straightforward manner. In the former case, the transfer matrix method must be 
modified. The results depend as well on which set of boundary conditions are chosen. 

The calculation of the Stark shift of the bound states in a single quantum well is described 
below. This example is chosen in order to illustrate the consequences of the form assumed by the 
applied field, whether it is uniform or localized. First, the solution for the unrestricted uniform field 
is outlined, assuming a uniform effective mass. Then it is explained how the solution must be 
modified when there is an effective mass jump at the well edges. Finally, the solution for a 
localized electric field is outlined. 
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2.B.6.a. Single Quantum Well In Unrestricted Uniform Electric Field 

2.B.6.a.i. Constant Effective Mass m e ff = mo 

Let us first consider the fate of bound states of a single well in a uniform electric field. 
(Refer to Figure 14) Recall that the bound states of a finite potential well are easily found by 
setting the element Mi i of the total transfer matrix to zero. This is the eigenvalue condition that 
follows from the requirement that the wave function be bounded on both sides of the well for 
energies below the top of the well. In that case the particle is localized in the well. Under a uniform 
electric field, however, true bound states no longer exist, because there is now the possibility that 
the particle can lower its energy by tunneling through the side of the well. There will still be 
energies, however, where the residence time of the particle in the well is long. These are close to 
the zero-field bound state energies, but are shifted by the field. The goal of this calculation is to 
find the quasi-bound states and their functional dependence on the applied field. 

The criterion for quasi-bound states is not that the wave function must be bounded on 
either side of the well, but that its amplitude be maximized in the region of the well. The wave 
function on the downstream side of the well is composed of an incoming and outgoing part: quasi- 
bound states exist when the ratio of the incoming to the outgoing parts is maximized. 

Refer to Figure 14 for the geometry of the well in an electric field. The Schroedinger 
equation in the regions in and around the well is: 

rJtL - e F z 4 > - Eifr |z| > a (69a) 

2 m eff dz 2 P 

^4 - <V C + C F z) * = E h I z | < a (69b) 

2 m eff dz 2 

Inside the well, the slanted potential can still be broken up into intervals of constant potential 
where the wave function is the sum of forward- and backward-propagating plane waves: 

i k: Z — ik;Z 

= Aj0 ' + Bje 1 

The wave-vector in each interval, kj, is easily calculated knowing the field and the coordinate at the 
edge of the interval: 

_ |70 ) 

1 ' * 
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The transfer matrix accumulated between the first and last intervals relates the coefficients of the 
wave functions: 


B 

c 


M„ • M 1 


• M 2 . . . M f 



(71) 


This matrix can now be used to relate the wave functions on either side of the well. 

Outside the well, where |z| > a, because the field is uniform, the Schroedinger equation is 
most easily solved by a transformation of coordinate. A new coordinate, x, is defined which is a 
dimensionless length minus a dimensionless energy, or: 


X 



2m eff E 

12 m eff e F ap^ )2/3 


(72) 


Inserting (-x) into equation 69a results in the Airy equation (31): 


d 2 

9 — x\p “ 0 (73) 

d x^ 


This has the general solution: 


i// = C Ai( -x) + D Bi (-x) (74) 

Figures 15 and 16 show the functions Ai(-x) and Bi(-x). Inspection of equation 73 for x shows that 
it is very large at small values of the field F. When this is the case, asymptotic forms of the Airy 
functions can be used (32): 


Ai (x) 

~ — j-j— 1/2 2 — 1/4 
2 

2 ( “ 1)k C k 
0 

| arg z | 

< II 

(75a) 



OO 




Bi (x) 

=* n ~ 1/2 z- 1/4 

2 C k f- k 

| arg z | 

<-!]- 

•x 

(75b) 


where: 

y = -L z 2/3 ( 76 ) 

3 
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The requirement of a bounded wave function on the left side of the well means the wave function 
there consists of Ai(-x) alone: 


\p — = a Ai ( — x) 


(77) 


On the right hand side of the well, the wave function is composed of both Ai(-x) and Bi(-x): 


i P + = E Ai (- x) + F Bi ( -x) 


Referring again to Figures 15 and 16, shows that the component to be maximized to the 
right of the well for a quasi-bound state is Bi(-x). If Bi(-x) is maximized at the expense of Ai(-x), the 
amplitude of the wave function is peaked near the well. The problem of finding the quasi-bound 
states is a matter of searching for the values of the energy E at a given field strength which 
maximize the ratio F/E. 

The ratio of the leftmost wave function to its derivative, at z = -a, gives the value of the 
wave vector there, ko: 


Ai ( -x) 


Ai' 


-x) 


z = — a 


(78) 


The transfer matrix Mq can now be initialized at z = - a: 



Now, matching the wave function and its derivative in the last interval of the well to that on 
the right at z = a allows the ratio of the incoming and outgoing waves, F/E, to be obtained in terms 
of the applied field: 


Pe ik f a + 


Qe 


— i kj a 


E Ai ( -x) + F Bi ( -x) 



(80a) 


Pj kf e 


kfa 


— Q: kz e 


— ik^a 


2fr>ef f f . E Ail— x) + FBi(-x) 


ti 2 


(80b) 
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Solving for the ratio F/E between the two equations above gives: 


F/E = 


Ai* (-x 2 ) 


l(-Je ik ' a _ ikf r ik f 8 ) - ^e ik » a +e - ik f a )] 


Bi- (-x 2 ) 

Ai' (— x 2 ) A 


_J /P J k f a + _~ ik f a \ Bi( ~ x 2 > _J_/_Z e ik,a +i 

i(— x 2 )\Q / Ai (— x 2 ) Aj' (— x 2 >7\ Q 


k f e 


- ikfi 


( 81 ) 


where: 



P/Q is determined by the total transfer matrix T, the coefficient of the leftmost wave function, A, 
and the wavevector at the left edge of the well, ko: 



The expression above can be put in terms of these parameters by matching at z = -a. The wave 
vector just to the right of the left edge of the well is: 


kj = V 2 m eff (E+V o-| F aP Ca l> 


( 83 ) 


and the wave vector just to the left of the right edge, at z = a, is: 


Inserting kj and kf into equation 81 and utilizing equation 78 for ko in the expression for P/Q 
(equation 82), will result in an expression for F/E explicitly in terms of the applied field Fap, incident 
energy E, well geometry, and total transfer matrix T. 

It is straightforward to write a computer program to calculate F/E from in equation 81 . The 
ratio F/E can be plotted as a function of the energy for varying field strengths, and the quasi- 
bound state energy levels appear as the maxima in this ratio. 

2.B.6.a.ll. Non-Constant Effective Mass: Jump In m«ff at the Well Edges 

In quantum wells formed of semiconductor heterojunctlons, the effective mass can 
change sharply at the edges of the well. When this happens, the modified connection rules must 
be used when matching the plane waves to the Airy functions at the well edges. If the effective 
masses inside and outside the well (m 2 and mi , respectively) are retained a new expression for 
P/Q results which differs from the one given in equation 82 only by a factor of (nrt 2 /mi) 1/2 : 



F/E is now found by inserting equation 85 into equation 81 : 



( 86 ) 


It should be noted that the effective mass difference also makes its presence known in 
the values of ko, kj, and kf. An effective mass jump at the edges of the well affects the number and 
energies of the calculated bound states in the well. 

2.B.6.b. Single Quantum Well: Localized Electric Field 
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Figure 17 shows a single quantum well in a localized electric field. The field results in a 
linear potential restricted to the areas inside and just on either side of the well. Outside this region, 
the field is zero and the potential is constant. This configuration represents real physical situations 
better than the unrestricted uniform field picture: in real semiconductor quantum wells and 
superlattices, the structure is biased by means of cladding layers in ohmic contact with metals at 
constant potential. Outside the cladding, the field is ideally zero. 

This is physically a different situation from the unrestricted uniform field case. There, as 
pointed out, no true bound states exist. When the field is localized around the well, however, 
bound and quasi-bound states may exist simultaneously. Whether or not they exist depends on 
the geometry of the well, the extent of the applied field, and the field strength. Figure 18 shows 
how either bound or quasi-bound states, or both, may arise. An obvious condition on the applied 
field, if true bound states are to exist, is that: 


v„ 

F < -2 (87) 

a 

Reid strengths higher than this allow only quasi-bound states. 

The bound and quasi-bound states of the well, and their Stark shift, are easily found by 
means of transfer matrices. Here the environs of the well are surrounded by constant potential 
regions. The slanted regions around the well are broken into intervals (see Figure 19) over which, 
once again, the solutions to the Schroedinger equation are the plane waves. Then, as outlined in 
section 2.B.5, for arbitrary potentials, the overall transfer matrix is found by multiplication. The 
criterion for a bound state is that Mu be zero. For a quasi-bound state, Mu need only be 
minimized. To calculate the Stark shift of a level, Mi i is calculated as a function of the applied field: 
the shift of the minimum under the field is the Stark shift of the level in question. 

2.B.6.C. Multiple Step Barriers: Constant, Localized Electric Field 

A solution for the transmission resonances of an MSB with N barriers, under zero external 
field, was outlined in a previous section. The solution was straightforward mainly because of the 
simple form of the zero-field Hamiltonian, and because of the right-left symmetry of the MSB. 
When a constant electric field is applied, that solution is no longer valid. The transfer matrix 
method is easily applied to calculating the Stark shift of the MSB’s resonances, and hence the 
tunneling current-voltage curve, when a constant, localized electric field is applied. Figure 20 
shows an N-barrier MSB, of overall length I, and barrier height Vq, under an applied field F ap . The 
barrier widths are a, and the well widths, w. In each barrier and well region, the Schroedinger 
equation is given by: 
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( 88 ) 




+ (v„ 


e F ap z) * = E * 


(barrier) 


-h 2 d 2 

2m w dz 2 


e F an z i// » E \p (well) 

ap 


( 89 ) 


Outside the MSB, the applied field is zero. The barrier and well regions of the I th period are 
broken up into nb and n w intervals of constant potential, respectively. Within each interval, just as 
for the isolated well, the solutions to the Schroedinger equation are the backward- and forward- 
propagating plane waves: 




= A 


I ,i,b,w 


J k l.i.b.w z 


+ B 


l,i,b,w 


-ik l i b w z 


( 90 ) 


The potential drops at the steps are all the same size because the field is constant inside the 
MSB: 

e F__a a F_ n w 

AV b = 2E- Ola) AV b = — ( 91b ) 

n b n w 


The transfer matrix at each step in the I th barrier or well is: 



and mb and m w are the effective masses in the barriers. and wells There are (nb,w - 1) of these 
matrices in each region. 

Transmission through the large steps at the right and left edges of the barriers is also 
handled by transfer matrices. The matrix through the left edge of the I 1 * 1 barrier is: 
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( 94 ) 



where: 


z L | = (I — i) (a + w) 


(95) 


and: 


with: 


‘l,i 


_\/2™h 


e-( v l.i +v q) 


(96a) 


- _ V 2 m w (E-V u ) 


L.l 


(96b) 


V |_ | = e F a p (I — 1) (a + w) (97) 

The transfer matrix through the right edge of the I th barrier is: 



where: 

Z R I = I (* + w ) — w (99) 


and: 


with: 


+ _V 2 m w (E-(V R 1 + V 0 )) 

RJ ' 1, 


(100a) 


- _V 2m w (E " V R.|> (100b) 

R 'l t. 


V R,I = eF ap (l (a + w)-w) 


( 101 ) 
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The overall transfer matrix is found, as usual, by multiplication. Then the transmission 
coefficient can be obtained as a function of both energy and applied field. This gives the Stark 
shift of the resonances. 

Tunneling current-voltage curves are obtained using equation 20: 

J(V,_,E) = — /— \ T (v ,E) 
kj \m j / a P 

where Vgp = -e F ap I, and E is treated as a parameter. In computing the J-V curve for an AIGaAs 
MSB, E can be taken as the Fermi energy in the n-cladding layer, typically a few thousandths of an 
eV. The effective masses mt and mj, are those of the n- and p-cladding layers. The wavevectors in 
these regions are kt and kj from above: 

v/2m: E -y/Zm (E — V ) (ln ~. j 

kj = V ^ (102a) k t = ^ ^ 5E_ (102b) 


It is easy to write a computer program which varies F ap (or V ap ), and for a given MSB 
geometry, calculates the J-V curve. Comparison of the J-V curve with the T-E plots showing the 
Stark shift of the resonances, clarifies the role that resonant tunneling plays in conduction 
through these structures. This is discussed in detail in a subsequent section. 

In AIGaAs MSB’s, where the effective mass varies between the wells and barriers, the 
prefactors of the transfer matrices at barrier edges cancel when the matrices are multiplied. The 
effective mass difference does survive, however, in the wavevectors, and should not be 
neglected. The effective masses in both cladding layers may differ as well, and their ratio occurs in 
equation 20, which gives the tunneling current. The magnitude of the tunneling current, the 
location of the resonances, and hence the shape of the J-V curve, are all affected strongly by the 
effective masses used throughout, and care should be taken to use the appropriate values. 

2.B.6.d. Sawtooth Barriers: Zero and Non-Zero External Electric Field 


Tunneling calculations for sawtooth barriers of the kind shown in Figure 21, under zero 
external field, require solution of a Schroedinger equation for the which the Hamiltonian is: 


H 




(103) 
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I 1 H Hill i 


where: 


e F eff = <V») ,104) 

is the slope of the barrier. This is effectively the same as the Hamiltonian for a step barrier or well 
under a constant, external electric field. The Schroedinger equation for the zero-field sawtooth 
barrier is: 


l > 2 

2 m eff dz 2 


+ eF eff zi It = 


(105) 


If an external, constant electric field is applied as well, another term, (e F ap z), for the 
external field, is added to the Hamiltonian and the Schroedinger equation becomes: 


ti 2 d 2 0/ 
2m eff dz 2 



(106) 


If the effective mass is not a function of z, a direct solution in terms of the Airy functions is possible 
in the same way it is carried out for a quantum well in a uniform field. Actual devices, however, are 
made from graded band-gap AI(x)Ga(1-x)As, and the effective mass is indeed dependent on z. By 
using the transfer matrix method, this complication of the Schroedinger equation is circumvented 
by approximating not just the potential, but the effective mass gradient as well, by a series of small 
steps of constant value. 

Suppose there is an N-tooth sawtooth barrier for which the T-E and J-V curves are 
sought. Each barrier is divided into n intervals of width Az = b/n, and for each of which the 
potential drop is given by: 



n 


The solution to the Schroedinger equation in each interval is, as usual: 


where: 


*lj (z) 


ik , :Z 

a *•» 


-ik, . z 

B u' 


(107) 


(108) 
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( 109 ) 


\/2m ef f (j Az) (E - (V 0 /b + e F gD )(l b + (j - 1) Az)) 

ti 


The transfer matrix at each step is: 

kj jj.i \ il^i ;-n “ ^1 z i 

l,i 


.-.if: 

l >' 2 N n 


m. = = — ha^L 


,, K yti\." K i.j*i- K i.i' z i.i (, k '.KiV~ l(k '.i^ + k '.i lz u 

k, 


k u 


i_HyiiV <k u + i + k u ,z u /i + ^yn\ e - ,lK u+i- K i' z j 


-i(k- -k:)Zi 


C U 




( 110 ) 


where: 


Mj +1 


_V 2m eff (z i + Az) ( E ” (V o /h + e F ap ) z l,j ) 


( 111 ) 


and: 


z, j = (lb + jAz ) 


(112) 


The transfer matrix through the right-most edge of the I th sawtooth handles a potential 
drop of Vo, the peak barrier height: 



I ~ R i 


where: 


i - -L 



(113) 


* _ -/Tm-IE- 


% 


(114a) 


■J 2m'j"(E — VJ ) 

" t. 


(114b) 


and: 


V T = < V o + eF ap z l> 


(115a) 


Z, = lb 


V | = e F ap z l 


(115b) 

(116) 


m , = m e ff (l b) “ m max 


(117a) 


m f = m eff (0) = m min 


(117b) 
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The transmission coefficient and tunneling current are computed from the element Mi i of 
the overall transfer matrix exactly as for the MSB described in the previous section. Assuming that 
the effective mass is graded in some manner over the barrier width, when the multiplication of all 
the transfer matrices is carried out, the prefactors containing the ratio of the effective masses to 
the right and left of each step, do not cancel as they did for the MSB with alternating layers. This is 
a consequence of the asymmetry of the sawtooth barrier. In AIGaAs sawtooth barriers, the 
triangular shape is accomplished by linear aluminum concentration grading. The attendant 
gradation of the effective mass must be accounted for by considering the functional dependence 
of m e ff(z ) at each step when the transfer matrices are computed. This is because the effective 
mass has a strong effect on the number and energy levels of the resonances supported by the 
structure, just as in the MSB and quantum well. 

The AIGaAs multiple sawtooth barrier which finds application as an avalanche photo-diode 
is operated under high reverse bias. This means that the applied field is of opposite sign to the 
effective field caused by the band-gap grading in each layer. As the applied field is raised, at some 
value it will exceed the band-gap grading and the sawtooth will become a staircase structure. 
(Refer to Figure 21) The value of the applied field where this occurs is: 


F 


staircase 



( 118 ) 


At this point there are no longer any barriers through which to tunnel. It is expected that the 
transition from sawtooth to staircase should be reflected by a change in the character of the J-V 
characteristic. In this respect a multiple sawtooth barrier is much different than a multiple step 
barrier. 


2.C. AI(x)Ga(1-x)As Heterostructures 
2.C.I. Properties of AI(x)Ga(1-x)As 

Superlattices of quantum wells and barriers, in real semiconductor materials such as 
AI(x)Ga(1-x)As, are made possible by the conduction and valence band discontinuities occurring 
at interfaces of materials with different band structures. In AI(x)Ga(1-x)As, the energy band 
structure depends on the aluminum content, x. By combining layers of differing Al content, it is 
possible to adjust the band-gap on either side of the interface and in that way to tailor the 
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conduction and valence band discontinuities. These discontinuities are what form the barriers to 
electron and hole transport through the superlattice. In tunneling calculations for AI(x)Ga(1-x)As 
superlattices, it is thus necessary to know the relationships between composition, energy band 
structure, and band alignment at heterojunctions. 

Carrier transport through a superlattice is also dependent on the carrier effective mass at 
every point in the structure. The effective mass in semiconductors is generally a function of the 
energy band structure, and in AI(x)Ga(1-x)As it is therefore dependent on x. Tunneling 
calculations in AI(x)Ga(1-x)As thus also require some knowledge of the relationship between 
composition, band structure and effective mass. 

In this study, only tunneling of electrons between direct conduction bands is considered. 
The information used in these calculations is: the direct conduction band alignment at GaAs- 
AI(x)Ga(1-x)As heterojunctions, and the electron effective mass in the direct conduction band 
minimum, as a function of x. The effect of non-parabolicity of the bands on the effective mass of 
electrons at high energy, for instance, when the electron encounters a step in a staircase 
superlattice, is neglected. The effective mass is assumed scalar. 

In the sections below, there is a brief discussion of the energy band structures of GaAs, 
AlAs, and AI(x)Ga(1-x)As, covering only the properties needed to perform tunneling calculations 
for AI(x)Ga(1-x)As superlattices in the effective mass approximation. 

2.C.1.a. GaAs Energy Band Structure and Effective Masses 

Pure GaAs is a direct-gap material: its conduction and valence band extrema occur at the 
Brillouin zone center. Figure 22 shows the <1 00> and <1 1 1 > energy band diagram for GaAs (33). 
The two lowest-energy conduction bands are shown: ri c , and Xi c , the direct and indirect bands, 
respectively. Also shown are the three highest-energy valence bands: Ti v and T 2 v. heavy and 
light-hole bands, respectively, and a third, T 3 V , split off by 0.36 eV from the other two by spin-orbit 
coupling. 

The direct band-gap, at the zone center, is 1.44 eV. The indirect X-| C minimum is 0.38 eV 
above the direct minimum. In thermal equilibrium, free carriers are found only in the lowest 
conduction band minima, although scattering of electrons by strong electric fields can excite them 
into X-| C . (29) This effect has been neglected in this study. Density-of-states effective masses for 
electrons and holes in GaAs have been determined (18), and are used throughout the present 
study: this is quite appropriate for the Ti c band in GaAs, for which the constant E(k) surface near 

the band minimum is spherical, and hence, for which the effective mass is actually scalar. (For the 
hole bands r 1v and T 2 v. this is a less satisfactory approximation, since these bands are non- 
parabolic and anisotropic.) The electron effective mass mp at Tic is 0.0636 ±0.002 (at 290 K); for 
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electrons in Xi c , mx is estimated as 0.39. The heavy-hole effective mass in Ti v is approximately 
0 . 68 . 


2.C.1.b. AlAs Energy Band Structure and Effective Masses 

Figure 23b (34) shows the <100> energy band structure of pure AlAs: inspection of this 
diagram shows that AlAs is an indirect-gap material, unlike GaAs. As Figure 23 shows, however, 
the symmetry of the conduction and valence bands is similar in the two materials. The indirect 
band-gap in AlAs is 2.2 eV; the band gap at the zone center is 2.9 eV. The energy separation of 
Tic and Xic is 0.877 eV. As in GaAs, there are two degenerate valence bands with maxima at the 
zone center. 

The density-of-states effective masses for AlAs have been determined as follows: in Xi c , 
mx is 0.37, in Tic mp is reported as 0.128 to 0.15, depending on the method of calculation. (18) 
The heavy-hole mass in Ti v is estimated as 0.85. 

2.C.I.C. AI(x)Ga(1-x)As Energy Band Structure and Effective Masses as a 

Function of x 

In the mixed crystal AI(x)Ga(1-x)As, the band structure and effective masses depend 
strongly on the aluminum concentration x. Since the pure compounds AlAs and GaAs are indirect 
and direct gap materials, respectively, there is naturally an aluminum concentration at which the 
AI(x)Ga(1-x)As changes from one type of material to the other. 

Figure 24 shows how the energy gap in AI(x)Ga(1-x)As depends on x at room 
temperature (19). The band-gap at Ti c increases with x, from 1.43 eV in pure GaAs, to 2.9 eV in 
pure AlAs. The experimentally measured band-gap at Ti c is (18): 

E o = 1.424 + 1.266 X + 0.26 x 2 (119) 

The quadratic term is relatively small compared to the linear term in x, hence Eo is roughly linearly 
dependent on x. The gap at Xi c increases less rapidly with x, and intersects the Ti c curve at a 
band-gap of 1 .96 eV, when x is approximately 0.4. The band-edge separation between X 1c and 
ri c has been reported as (18): 


A = 0.380 - 0.892 x - 0.365 x 2 


( 120 ) 
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For x greater than 0.4, AI(x)Ga(1-x)As is an indirect-gap material. In fact, although the X 1c - 
r 1c crossover occurs at x = 0.4, Figure 25 (18) shows that already half of the conduction electrons 
are in the Xi c minimum when x is only 0.3. This is because, even though the Xi c minimum is still at 
higher energy than that of Fi c , the density of states is a factor of 45 greater at Xi c (18). This 
means that for GaAs-AI(x)Ga(1-x)As heterojunctions where x is greater than 0.4, direct-indirect 
band tunneling plays a dominant role in carrier transport across the interface. Since the 
calculations presented here address only direct-direct band tunneling, treatment of GaAs- 
AI(x)Ga(1-x)As heterojunctions is restricted to compositions with x < 0.4. 

The density-of-states effective masses for AI(x)Ga(1-x)As, obtained by linear interpolation 
with x (18), gives the following relation for mx, the effective mass of electrons in Xi c : 


m x = 0.39 (1 - x) + 0.37 x 


( 121 ) 


For mp in Tic. the experimentally measured band-gap Eo was used to find the spin-orbit 
splitting of the mixed crystal; mp is estimated using this information as (18): 

mp = 0.0636 + 0.0552 x + 0.0092 x 2 (122) 

This is the relationship used throughout this work for the electron effective mass in AI(x)Ga(1- 
x)As: the effective mass mp ranges from 0.0636 in pure GaAs, to 0.0872 in AI(0.4)Ga(0.6)As. For 

sawtooth barriers in AI(x)Ga(1-x)As, the linear grading of Al results in approximately linear band- 
edge grading. Inspection of equation 119 for the band-gap as a function of x, however, shows 
that there is a slight bowing, indicated by the small quadratic term. Likewise, the linear aluminum 
grading results in a parabolic z-dependence of the effective mass m in the sawtooth barrier, 
reflected by the small quadratic term in equation above. The bowing of the band-gap was 
neglected, but the parabolic term in the effective mass was taken account of in these tunneling 
calculations for AIGaAs sawtooth superlattices. 

2.C.1.d. Band Alignments at GaAs-AI(x)Ga(1-x)As Heterojunctions 

The phenomenon giving rise to barriers to electron and hole transport through 
semiconductor superlattices, is the misalignment of conduction and valence band edges when 
materials with different band structures are brought into contact. The theoretical and experimental 
determination of these band misalignments in heterojunctions is difficult and, for GaAs-AI(x)Ga(l- 
x )As interfaces, is still not completely settled (13,35-40). In this study, the average values of the 
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best available measurements were used for the Ti c band-edge alignments. These are 
represented in Figure 26 (courtesy of R. C. Miller, ATT-Bell Labs) showing AEc as a function of x. 
Note the shift to the indirect minimum at x = 0.4, where AEc is 0.33 eV. Up to this point the ratio of 
the band-edge misalignments, AE C :AE V , is approximately 60:40. The 0.33 eV direct conduction 
band-edge misalignment at x = 0.4, is the largest attainable in GaAs-AI(x)Ga(1-x)As 
heterojunctions, before interband tunneling takes over. This may be a practical limitation on the 
efficiency of electron impact ionization in AIGaAs staircase avalanche photodiodes, since the 
probability of impact ionization at an interface is directly related to the magnitude of AE C there (3). 

2.C.2. Some Previous Tunneling Calculations In AIGaAs 

Heterostructures: Tunneling In Multiple Step Barriers, Stark Shift of 

Bound and Quasi-Bound States In Single Quantum Wells 

In the section below, some previous solutions to two important problems are outlined. 
Both problems are relevant to the problem of tunneling in graded band-gap superlattices. 

The first work to be discussed consists of tunneling calculations by Esaki and Tsu (11) for 
multiple step barriers formed by GaAs-AI(0.5)Ga(0.5)As heterojunctions. They obtained the 
transmission coefficient T(E), and the tunneling current J(V), for a small number of barriers. The 
approach they adopted is very similar to the transfer matrix method used in this study. Their results 
for step barriers are qualitatively somewhat similar to those expected for a sawtooth array of the 
same barrier height, composition, and roughly the same spacing. For these reasons, the results of 
Esaki and Tsu for MSB's were reproduced by the method used in this study in order to verify its 
accuracy, since there are no prior tunneling calculations for sawtooth superlattices with which to 
compare. 

The second work to be discussed is a calculation by Austin and Jaros (41), of the Stark 
shift of quasi-bound states in a single GaAs-AI(0.5)Ga(0.5)As quantum well, in an unrestricted 
uniform electric field. This work is reviewed primarily to show how the spatial extent of the applied 
field affects the calculated Stark shift. This is relevant to tunneling calculations for sawtooth 
superlattices, since practical devices will be operated by applying an electric field that is confined 
to the superlattice between two metal contacts. Assumption of an unrestricted uniform field will 
lead to the wrong functional dependence of the Stark shift. In the present work, the results of 
Austin and Jaros for a single quantum well in a uniform field, are reproduced by means of the 
hybrid Airy function-transfer matrix method detailed in section 2.B.6.a. above. Then, the same 
calculation is repeated for the well in a localized electric field. The results, which show the 
consequences of the geometry of the electric field for the calculated Stark shift, are presented in 
a later section. 
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2.C.2.a. Esaki - Tsu Tunneling Calculations for the Multiple Step Barrier 

Esaki and Tsu (11) performed tunneling calculations for conduction electrons in two, 
three, and five barrier GaAs-AI(0.5)Ga(0.5)As MSB's. The barrier heights were taken as 0.5 eV. 
(The data in Figure 26 giving aE c as a function of x, would have set the direct conduction band 
misalignment at 0.4 eV at x = 0.5, instead of 0.5 eV) The barriers were 2 nm wide and 5 nm apart. 
(See Figure 27) Working within the effective mass approximation, they do not cite the values used 
for the effective masses in the well and barrier layers, nor do they quote a source. It is likely that 
the values they used were close to those calculated by means of equation 1 22 for the r i c 
effective mass: 0.0636 in pure GaAs, and 0.0935 in AI(0.5)Ga(0.5)As. Thes e a re the values used 
in reproducing their results. In calculating the J-V curves for these MSB’s, Esaki and Tsu 
postulated a constant electric field distributed over the length I of the superlattice. They use a 
transfer matrix method which is outlined below. 

The total energy E for an electron in this one dimensional superlattice is the sum of 
longitudinal (in the direction of the superlattice) and transverse parts: 


E 


E C 



(123) 



i k R z 

= T^ t e (125b) 

R and T are the reflection and transmission coefficients for the structure. They are obtained from a 
product of the transfer matrices at each barrier: 
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where M p is obtained by matching y p and y' p at each interface: 



and where: 


_ 1 /2"'r.o <v p- EI 


(128) 


V p and mr, p are the potential and the effective mass in the p ,h layer of the superlattice. This 
transfer matrix gives the same results as the transfer matrix defined in section 2.B.6, however, the 
identification of the coefficients of the wave functions as in equation 125 above results 
immediately in the expressions below for the transmission and reflection coefficients: 


R 



(129a) 


T = Mii -M 


M 


’ 12 ‘ 
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(129b) 


This is the expression for T(E) that the authors use. 

In calculating the J-V curves for the MSB, they neglect the potential gradient within each 
layer, and take account only of the potential difference between the centers of each layer. (Refer 
to Figure 27) This is an approximation which, because of the small dimensions involved, is 
probably acceptable. 

The expression below is used to calculate the tunneling current: 

0° oo 

J=-£— f dkp f dk t ff(E) — f(E*)] T*T — £ (130) 

4rr 3 t, / / 3k *! 

Jo Jo 

where f(E) is the Fermi distribution function and: 

/ (Eo- E f )/kT. \-i 

f(E) = \e V + 1 / (131a) 
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(131b) 


f(E') 


= ( e (Eg-(E f -eV))/kT + ^ 


-1 


The Fermi level, Ef, is taken as 0.005 eV for GaAs with n = 10 17 /cm 3 . An integration of equation 
130 over the transverse direction using the Fermi functions in equation 131, gives: 


j = 


emp kT j 

I 

J n 


:f / 1 + P (E f- E e ,/kT 

T * T £n / 1 + e 


1 + e 


(E f — Eg — eV)/kT 


dE c 


(132) 


For small values of T, J is reduced to: 


j = 


e m T 


2t r 


2f, 3 


f f (Ef — Eg) T* 
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V > Ex 
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(133a) 


(133b) 


In this manner the contribution of all electrons within about ±kT of the Fermi level of the cladding 
layer, is integrated into the total tunneling current. 


2.C.2.b. Stark Shift of Quasi-Bound States In GaAs-AI(0.5)Ga(0.5)As 
Single Quantum Well 


Austin and Jaros present an exact numerical calculation of the Stark shift of the quasi- 
bound states in a single quantum well. The well they consider is one which they find supports a 
single bound state at zero field strength. The well is 3 nm wide, 0.4 eV deep, and is formed 
between two GaAs-AI(0.5)Ga(0.5)As heterojunctions. The authors work within the effective mass 
approximation, taking the effective mass of the conduction electrons to be 0.067 inside and 
outside the well, and neglecting the discontinuity at the edges. They find the Stark shift for both 
electrons and holes, but only the results for electrons are discussed here. 

The approach adopted is similar to the hybrid transfer matrix - Airy function solution 
outlined in a previous section. The Hamiltonian used is: 


h = 


~2 

JE_ 


2m 


eff 


- eF ap z 


+ V, 


(134) 
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(69a) 


for a uniform field Fap. The Schroedinger equation is then given by equations 69a and 69b: 

- — ~ eF ap z ^ = 1*1 >a 

2 mp dz^ 

_ _hL d2 }L - (V 0 + e F ap z) <// = Et|/ |z| < a (69b) 

2 mp dz 2 

where, as formerly, a is the well width and Vo its depth. The solutions to the Schroedinger 
equation in each region are: 

\p~ - aAi<— x — ) z < — a (135a) 


w = CAi (-x^ + DBi(-xJ 


| z| < a (135b) 


= E Ai (— x + ) + F Bi (— x + ) 


z > a (135c) 


where: 
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(136a) 


(136b) 


The authors match these wave functions and their first derivatives at both edges of the well, 
obtaining the coefficients of each of the components. They state that resonances, or quasi- 
bound states, are characterized by an abrupt increase of ji in the phase <|> in the equation: 


tan 



(137) 


This is another way of stating that the ratio of the downstream wave function components, Bi(-x ) : 
Ai(-x ), must be maximized for a quasi-bound state. By varying the applied field Fgp, calculating the 
ratio F/E, and searching for shifts of re, they obtain the quasi-bound state Stark shift. They find a 
Stark shift with quadratic dependence on the field, for field strengths between 100 and 500 
kV/cm. Their results agree with those previously obtained by means of variational calculations. 
(42,43) 

The authors neglect the effective mass difference inside and outside the well. Outside 
the well, using equation 92 and x = 0.5, m e ff is 0.0935, almost 50% higher than the value inside 
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the well. Neglecting this difference leads to an error in the calculated energy level of the 
resonances. This is discussed in greater detail in a later section, where their results are compared 
with those obtained in this study. 
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3. RESULTS AND DISCUSSION 


3.A. Scope 

3.A.I. Description of Work 

The aim of this work is to gain an understanding of the tunneling properties of electrons in 
graded band-gap AIGaAs sawtooth superlattices. All computations are carried out by means of the 
transfer matrices described in the theory section, working entirely within the effective mass 
approximation. T-E and J-V curves are here presented for these sawtooth structures. 

In order to put these results in the context of previous work on conventional step barrier 
superlattices, tunneling calculations for some other simple structures are also presented. Transfer 
matrix tunneling calculations for single and multiple step barriers under constant external electric 
fields are reported. The effects of the external field on the resonances supported by these 
structures are discussed, and J-V curves are presented for the MSB. In particular, results obtained 
here for AI(0.5)Ga(0.5)As-GaAs MSB's are compared to those obtained by Esaki and Tsu (ref) for 
identical structures. This is done to demonstrate the accuracy of the computational method used 
here. 

Bound and quasi-bound states of AI(0.5)Ga(0.5)As-GaAs single quantum wells are 
calculated, along with their Stark shift under constant electric fields. These calculations are done 
for the two configurations of uniform and localized field, in order to illustrate that, which of the 
above is chosen will determine whether the calculated Stark shift is quadratic or linear. These 
results are compared with those obtained by Austin and Jaros (ref) for an identical quantum well 
under uniform fields, again by way of validation of the present work. 

Finally, tunneling calculations for sawtooth barriers and superlattices are reported. The 
properties of sawtooth and step barriers are compared, and similarities between MSB's and 
sawtooth superlattices are pointed out. Throughout, particular emphasis is placed on the 
consequences of the effective mass variations in all the structures studied here, and need to use 
modified wavefunction connection rules at interfaces where the effective mass makes a 
discontinuous change. 

3. A. 2. Approach 

All computations are carried out by means of ihe matrix methods described in the 
preceding theory section. VAX Fortran computer programs computing the T-E and J-V curves for 
all step and sawtooth barriers and arrays, and the bound and quasi-bound states and Stark shifts 
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for single quantum wells, are found In Appendix A, along with descriptions of the operation of 
each program. All programs involving variations in the effective mass use the modified connection 
rules described in section 2.A. 

3.B. Single Step Barrier 

Tunneling calculations for single step barriers are done in the programs ZIGGURAT and 
ZIGGEFFMASS. The former assumes a constant effective mass equal everywhere to m 0 , while 
the latter models AI(x)Ga(1 -x)As-AI(y)Ga(1 -y)As step barriers, accepting the aluminum 
concentrations x and y as inputs, and calculating m e ff inside and outside the barrier from equation 
122. Both programs accept the zero-field barrier height Vo and width a, and the use of 
ZIGGEFFMASS requires choosing V 0 to agree with AE C given in Figure 26 as a function of x and y 
at the interfaces. The applied field F ap enters as a parameter: each run of the program then 
produces a file of T versus E, for a given value of F ap , which can be plotted. Both programs model 
the barrier in a localized field restricted to the barrier. 


3.B.I. Single Step Barrier: Zero external field 

Consider a single GaAs-AI(0.4)Ga(0.6)As step barrier 10 nm wide. The barrier height V 0 
must be 0.33 eV, according to Figure 26. The effective mass inside the barrier is 0.0872, while 
outside it is 0.0636. The T-E curve for this barrier is shown in Figure 28, where the energy ranges 
from 0.5 to 1 .5 times the barrier height. The first transmission maximum occurs at 0.372 eV, well 
above the barrier height: this curve is described by equation 28 for E > Vq, and by equation 31 for 
E < V 0 . 

To see the strong influence the effective mass can have on the T-E curve of step barriers, 
refer to Figure 29, which shows the curve for a barrier of identical height and width as the one 
above, but with constant effective mass equal to mo- The first maximum occurs at 0.334 eV, only 
slightly above the barrier height, and five more closely-spaced peaks follow. This curve is 
described by equations 29 and 32, above and below V 0 , respectively. Equation 30 gives the 

energies of the maxima in this case. 

Further differences between the two barriers are evident in Figure 30, where their T-E 
curves are overlapped. The effective masses mi and m 2 , and their relative difference, (mi - m 2 ), 
entering into equations 28 and 31, result in higher transmission minima throughout the range, in 
broadened resonances shifted to higher energies, and in lower-frequency transmission 
oscillations. The latter is the dominant effect, and is caused by the effective mass inside the 
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barrier, m 2 , appearing in the argument of the sine, being a great deal smaller than mo- The 
elevated minima however, reflect the effective mass difference (mi - m 2 ), which is negative and 
makes the denominator in equation 28 smaller. Figure 30b shows this: the T-E curve for a step 
barrier with constant m e ff = m 2 , is superposed on Figure 28. The two curves agree except in the 
absolute values of the transmission coefficient, which are highest when the mass difference (mi - 
m 2 ) is retained. In AIGaAs step barriers, this difference is always negative: hence the minima are 
always higher than those where the effective mass is taken as constant. 


3.B.2. Single Step Barrier: Non-Zero Field 

In considering the effects of a constant electric field on the transmission properties of 
single step barriers, it is easiest to start with the constant mass case (meff = mo), and then show 
how the effective mass further alters the picture when dealing with AIGaAs. Unlike the zero-field 
case, there is no analytical expression forT(E) for a step barrier under a constant field. 

T-E curves for the 10 nm, 0.33 eV high barrier were obtained, using the program 
ZIGGURAT, for applied field strengths F a p, of 2 and 5x1 O' 2 eV/nm. These are shown together 
with that of the zero-field barrier, in Figure 31. The field deforms the barrier to a thinner, triangular 
shape, and this facilitates tunneling at energies below Vq. In fact, resonance peaks occur at 
approximately 0.297 and 0.320 eV for F ap * 2x1 O' 2 eV/nm, and at 0.255 and 0.300 eV for F ap * 
5x1 0‘ 2 eV/nm. This is a major effect of the applied field, since resonance maxima below the top of 
the barrier are not possible when there is no external field, as equations 31 and 32 show. 

The broadening of the peaks evident in Figure 31, is a consequence of the tapered 
thickness of the barrier: for E greater than the lower edge of the barrier, there is no longer a single 
barrier width to serve as the resonance criterion: 

sin l <2 a = n IT (138a) 

Instead, there is a spread in the wavelength over which some kind of resonance can take 
place, because the width of the barrier varies. The field not only broadens the peaks, but causes 
them to be more widely spaced for the same reason. 

Another striking effect of the field is to lower the contrast between the resonance peaks 
and valleys. For any non-zero field, the peaks no longer reach 1, and the minima are raised with 
respect the zero-field barrier. The transmission curve asymptotically approaches 1 at infinite 
energy. As Figure 31 shows, increasing the field strength washes out the contrast of the peaks 
and valleys, and slows the approach to the maximum value of 1, while increasing the amount of 
tunneling below Vq. 
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When the 10 nm AI(0.6)Ga(0.4)As-GaAs step barrier is subjected to a constant electric 
field, effective mass variations further change the T-E curves. T-E curves generated by 
ZIGGEFFMASS for F ap = 2 and 5x1 O' 2 eV/nm are found in Figure 32, along with the zero-field 
curve. The increasing electric field has an even more pronounced effect on the transmission 
through the GaAs-AI(0.4)Ga(0.6)As step barrier, than it had when m e ff = mo throughout. The 
same trends are evident in much greater force: the peaks are broadened and shifted to lower 
energies, and their contrast is much reduced, relative to the zero-field case. 

Figure 33 compares the T-E curves of the constant mass and AI(0.4)Ga(0.6)As barriers. 
Figures 33a and 33b give the T-E curves for F ap = 2x1 0" 2 and 5x1 O' 2 eV/nm, respectively. Most 
of the structure in the T-E curves of the constant mass barrier is obliterated by the smaller effective 
mass of the AI(0.4)Ga(0.6)As barrier. Field-assisted tunneling is more probable through the 
AI(0.4)Ga(0.6)As barrier at all energies, at F ap = 5x1 O' 2 eV/nm. At F ap = 2x1 O' 2 eV/nm, this effect 
is less pronounced, but transmission minima for the AI(0.4)Ga(0.6)As are still higher than those of 
the constant mass barrier. 

Field-assisted tunneling calculations for AIGaAs step barriers should thus certainly use at 
least the correct average effective mass, and optimally, the effective masses appropriate to the 
barrier and adjacent layers. The former will result at least in the correct peak locations, although the 
absolute values of the transmission coefficient will be somewhat in error, away from the resonance 
maxima. Simply using the free electronic mass throughout will lead to gross error in the T-E curves 
obtained at all field strengths. 

3.C. Single Quantum Well: Bound States and Stark Shift 

In this section, the bound states of AIGaAs single quantum wells (SQW’s) and their Stark 
shift in a constant electric field, are obtained by means of transfer matrices. First, the zero-field 
solution for the bound states is given: the role of the effective mass in determining the energy 
eigenvalues of the SQW is described. The programs for the zero-field SQW bound states are 
WELL and ALGAASWELL. Both programs accept values for the well width a and depth Vq, and 
produce files of the total transfer matrix element Mi i as a function of the energy E below the top 
of the well, which is taken as zero. Minima in Mu indicate the bound states of the well. The 
program WELL assumes a constant effective mass m e ff = mo everywhere inside and outside the 
well. The program ALGAASWELL, however, models SQWs formed by AI(x)Ga(1-x)As- AI(y)Ga(1 - 
y)As heterojunctions, and, just as the program ZIGGEFFMASS does, requires the x and y values 
in order to calculate the effective masses inside and outside the well. Vq must be chosen in 
accordance with AE c (x,y) read from Figure 26. 
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The Stark shift of the energy eigenvalues of AIGaAs SQW's is obtained for constant 
electric fields localized to an area surrounding the well, and also for unrestricted uniform fields. 
This is done in order to show that, which of the two configurations is used, will determine whether 
the Stark shift is linearly or quadratically dependent on the field strength. 

The localized-field Stark shift is calculated in the program STARKMASS. This program 
models a 3nm wide AI(0.5)Ga(0.5)As well under the field as shown in Figure 18c. This is identical 
to the SQW studied by Austin and Jaros (41). The extent of the field is such that level of the 
plateau downstream of the well is at the same level as the lowest potential in the well. This makes 
sure that bound states in the well can tunnel out under a non-zero field. This allows a direct 
comparison with the results for the unrestricted uniform field, where all bound states under any 
field are susceptible to tunneling. In other respects, this program works identically to 
ALGAASWELL, except that the applied field Fap enters as a parameter. Each run of the program 
results in field-shifted minima in Mi i. Plotting the minimum Mu as a function of F ap then gives the 
Stark shift. 

The program AIRYWELL is used to calculate the Stark shift of a 3nm GaAs- 
AI(0.5)Ga(0.5)As SQW, under uniform electric fields. This program accepts the applied field as a 
parameter, and calculates the ratio (F/E) (of the inbound to outgoing wavefunction components 
downstream of the well), as a function of the energy E below the top of the well. Maxima of |F/E| 
indicate the bound state energy levels. Each run of the program gives these Stark-shifted levels 
at a given field. These levels are plotted as function of the field strength. 

3.C.I. Single Quantum Well: Zero-Field 

The program ALGAASWELL gives the spectrum for a GaAs-AI(0.5)Ga(0.5)As SQW, 3nm 
wide and according to Figure 26, 0.4eV deep. This is shown in Figure 34, where the total transfer 
matrix element Mu is plotted as a function of of E. This well supports a single bound state at - 
0.197 eV. This is in agreement with the results of others (ref) who first determined that these well 
dimensions, along with an effective mass of 0.0636, result in a single level. There is also rough 
agreement with an approximate calculation of the number and level of bound states, using 
equations 66b and 66c. These predict a single level at -0.211 eV. 

The discrepancy between these values is a consequence of neglecting the effective 
mass difference between the well and the adjacent layers. The program ALGAASWELL uses the 
appropriate values of 0.0636 outside, and 0.0935 inside the well. If instead, an effective mass of 
0.0636 is used throughout, then a single bound state at -0.217 eV is obtained. This is in excellent 
agreement with an easily derived equation for the energy of the ground state of a finite quantum 
well (44): 
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Figure 35 shows the Mu versus E curves overlapped for these two treatments of this 
well: the bound state level is strongly affected by how the variation of m e ff at the well edges is 
handled. This illustrates the risk involved in applying standard formulae to semiconductor 
heterostructures. Most, if not all, of these take no account of the effective mass variations at 
interfaces, and this will impair the accuracy of the results obtained. 

To illustrate further how profoundly the effective mass influences the bound states of the 
well, calculations for an SQW 3nm wide and 0.4 eV deep, but with the mass constant and equal to 
mg, were repeated in the program WELL. The results are shown in Figure 36. Instead of a single 
level there are four, and this is in agreement with the equation below giving the number of bound 
states of a finite quantum well (44): 


N >^ 2meffN ^ > n — 1 ... _ (139) 

7rti 

Inspection of this equation shows that the number of levels that a well can trap is 
proportional to the square root of the mass of the particle in the well. Since mo is approximately 
sixteen times as great as m e ff in GaAs, there are four times as many trapped states in this well. 
Again it is apparent that the appropriate effective masses must be used in calculations intented to 
model real semiconductor heterostructures, potential wells in this instance. 


3.C.2. Single Quantum Well: Stark Effect 

When an electric field is applied to a single quantum well such as the one described 
above, the bound states shift to lower energies inside the well. The magnitude and field-strength 
dependence of this energy shift depend on whether the field is localized or unrestricted. 

3.C.2.a. Single Quantum Well: Stark Effect Under Localized Field 

First, the Stark shift for the 3nm AI(0.5)Ga(0.5)As single quantum well in a localized field is 
calculated. In order to make a direct comparison with the work of Austin and Jaros and others (41 - 
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43) the effective mass is taken as 0.0636 everywhere, neglecting the discontinuity at the well 
edges. It was shown above that this results in an error in the determination of the bound state 
energy. It has no effect on the functional dependence of the Stark shift, however. In the program 
STARKMASS the total transfer matrix element Mi i is computed as a function of the energy and 
field strength. These data are plotted in Figure 37; the field strengths range between 0 and 5x10* 
2 eV/nm. The bound states, just as for zero field, are the minima in Mi It is easy to see the shift of 
the bound state energies to deeper in the well as the field increases. 

When these shifts are referred to the zero-field bound state energy of -0.21 7eV with the 
differences plotted in Figure 38 as a function of the field strength, it is clear that the Stark shift is 
linear. This effect is either explicitly described in the work of others, for instance, Schwartz in his 
study of AIGaAs MSB's (17), or is implied in the results of still others. Esaki and Tsu (11) for 
example, report J-V characteristics for MSB's, which, interpreted in the light of T-E curves for the 
same structure, indicate a linear Stark shift of the resonances. (This is described in detail in the 
following section). Here, the Stark shift of an isolated well is calculated. The common feature of 
the isolated well and the MSB's studied in the works referred to above, is that the field is localized 
around the structure in question. 

3.C.2.b. Single Quantum Well: Stark Effect Under Uniform Field 

Now the same isolated AI(0.5)Ga(0.5)As well in a uniform, unrestricted electric field is 
considered. The program AIRYWELL computes the ratio F/E of the well as a function of energy 
and applied field strength. Figure 39 shows F/E versus E for several field strengths between 
1x1 0‘ 2 and 5xl0' 2 eV/nm. The peaks in F/E, analogously to the minima in Mu, indicate the 
bound states. The shift in the location of the bound state is again to values deeper in the well, but 
this time the shift is obviously not linear. 

Figure 40 shows the bound state energy levels plotted as a function of the field strength. 
A polynomial least-squares fit to the data results in a quadratic Stark shift: 

E,„ = -221.85 - 0.9217 F* (140) 

W «K 


The energy shift with field strength is plotted in Figure 41 along with the same data for an 
identical well reported by Austin and Jaros. Fairly good agreement exists between the two sets of 
data. A larger shift is obtained by the present means, however, and the source of the discrepancy 
is uncertain. 
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The Stark shifts obtained from the localized and unrestricted fields appfied to the well are 
plotted together in Figure 42. The linear shift for the localized field is much larger at a given field 
strength than the quadratic shift. In fact, the magnitude of the linear Stark shift is calculated 
assuming that the geometric origin is centered at the left hand edge of the well. Comparing the 
linear Stark shift with the product of the field strength and the well half-width, strews that the Stark 
shift is entirely made up of the lowering of the potential in the center of the well, under the field. 
This is also observed when the Stark shift of the resonances of MSB's subjected to local fields is 
calculated. The linear shift under localized fields appears actually to be just the shift in the base of 
the well. For MSB's and isolated quantum wells alike, the virtual state is then still the result of 
resonance of the wave off the sides of the well, only the wavevector is referred to the bottom of 
the well, which is now shifted lower by the field. 

The Stark shift under unrestricted fields is in agreement with the results of others (41 -43). 
It is reported, however, that under stronger electric fields, (i.e. greater than about 1x1 O' 1 eV/nm) it 
departs from its quadratic dependence on the field. This regime is not explored in the present 
study. — 

In any case, it is clear that calculations of the Stark shift of the resonances in actual 
semiconductor heterostructures should probably be made by assuming that the electric field will 
be localized to the structure. The consequences of assuming an unrestricted field are that the 
wrong functional dependence and magnitude will result. For this reason, the field-assisted 
tunneling calculations for AIGaAs sawtooth superlattices, to be described below, are all performed 
assuming localized fields. 


3.D. Multiple Step Barriers 

Tunneling calculations for multiple step barriers are presented in this section. The 
transmission coefficient is computed as a function of incident energy, just as for the single step 
barrier, under zero external electric field. In addition, the tunneling current as a function of applied 
voltage is calculated. 

The tunneling calculations for AIGaAs MSB's are performed in the programs ALMSB and 
JVALMSB. ALMSB computes the transmission coefficient for an MSB of an arbitrary number, N, 
of periods, where a period is composed of a barrier and the well to the right of it. The applied 
electic field may be varied as a parameter, so that T-E plots parametric in Fgp may be generated for 
a given MSB. The program is similar to ZIGGEFFMASS, used for single AIGaAs step barriers, 
except for the fact that the total transfer matrix is accumulated through several barriers instead of 
one. 
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JVALMSB computes the tunneling current as a function of the applied voltage, creating 
data files from which J-V curves are made. It is a derivative program of ALMSB, differing only in the 
use of a single incident energy of 0.005 eV, roughly equal to the Fermi energy in bulk GaAs. The 
applied voltage is in this case the independent variable. 

Both of these programs have analogues which compute the transmission coefficients and 
tunneling currents through ideal structures where the effective mass is constant and equal to mo. 
These are the programs MSB and JVMSB, respectively. Except for using a single, constant 
effective mass throughout, they are identical to ALMSB and JVALMSB. 

In the text below, first the T-E and J-V data are reported for 2, 3, and 5 barrier GaAs- 
AI(0.5)Ga(0.5)As MSB's. These are the same structures studied by Esaki and Tsu (11): the results 
of the present study are described and compared with theirs. The correlation of the transmission 
resonances with features in the J-V characteristic is discussed. 

Then, the same calculations are repeated for GaAs-AI(0.4)Ga(0.6)As MSB's of the same 
dimensions and number of periods. This is done to provide a basis for comparison with the GaAs- 
AI(0.4)Ga(0.6)As sawtooth tunneling calculations in the final section. 

Results of tunneling calculations for ideal, m e ff = mo, MSB's are presented as well, in 
order to show once again the strong influence the effective mass has upon the resonances of a 
quantum heterostructure. 


3.D.I. GaAs-AI(0.5)Ga(0.5)As MSB's: Comparison with Results of Esaki 
and Tsu 


3.D.1.a. GaAs-AI(0.5)Ga(0.5)As MSB: Transmission Coefficient versus 

Energy 

The MSB's for which Esaki and Tsu did tunneling calculations consisted of 2, 3, and 5 
step barriers 2 nm wide and 5nm apart, and 0.5 eV high. These are the dimensions used in this 
study .The barrier height of 0.5 eV is not an accurate representation of Ec for this composition. A 
value of 0.4 eV, which was used in the bound state calculations for the quantum well, is more 
appropriate, according to Figure 26. The barrier height of 0.5 eV is used in the tunneling 
calculations presented here, however, for purposes of comparison. 

The transmission coefficients calculated in the program MSB, are plotted as a function of 
energy in Figure 43a, b and c, for 2, 3 and 5 barriers respectively. Below the barrier height in each 
case, there are two sets of resonances. The two-barrier MSB has a single resonance at 0.12 eV 
and another, much broader, one at 0.43 eV. 


55 


As pointed out in section 2.B.3, the occurance of resonances below the barrier height 
(under zero field) is impossible in tunneling through a single step barrier. When two or more 
barriers are involved, resonance of the wave with the wells between them is possible at energies 
much below their height, and that is what is giving rise to the resonances seen in Figure 43. 
Furthermore, each resonance peak reaches 1 . As demonstrated in section 2.B.3, this is 
characteristic of superlattices with even spatial symmetry, and these MSB's have this symmetry. 
When it is broken by applying an electric field, the peaks no longer reach 1. 

The three- and five-barrier MSB's have resonances which are split by coupling between 
the two and four wells, respectively, formed by the barriers. The T-E curves for the two- and five- 
barrier MSB's are shown on the same scale in Figure 44. There are four resonances distributed 
around the each single peak belonging to the two-barrier structure; two at higher, and two at lower 
energy. The off-resonance transmission coefficient is several orders of magnitude lower at most 
energies for the five-barrier MSB due to the greater difficulty of tunneling through three more 
barriers. 

The T-E curves calculated by Esaki and Tsu for each of the structures, are found in Figure 
45. These are in good qualitative agreement with those shown Figure 43, however, the 
resonances are now centered at 0.08 eV and at 0.32 ev, instead of 0.12 eV and 0.43 eV. The 
source of this discrepancy is unknown, however, the location of the resonances is very sensitive 
to the values of the effective mass that are used in these calculations. It is likely that the authors 
used somewhat different values than those used in this work; since they do not cite their values, 
this is unresolved. Agreement is good enough to conclude, however, that the method used in 
this study is vald. 

3.D.1.b. Tunneling Current versus Applied Voltage 

The tunneling current as a function of the applied voltage is calculated for these GaAs- 
AI(0.5)Ga(0.5)As MSB's in the program JVALMSB. The incident energy is 0.005 eV, that is, the 
Fermi level in n-doped GaAs where n = 10 17 /cm 3 . The applied voltage ranges between zero and 
2 volts. This range is sufficient to show how the resonances of these MSB’s influence their J-V 
characteristics. 

In Figure 46 are found the J-V curves for the two- and three-barrier MSB’s. The curve for 
the two-barrier MSB is particularly simple because there is only one well in the structure. There is a 
current peak at approximately 0.23 V, and another broader one at 0.86 V. On the high-voltage 
sides of these peaks are regions of negative differential resistance. These features are 
consequences of the resonances supported by the MSB, and are manifestations of resonant 
tunneling. There is a current peak at 0.23 V because half of this voltage is dropped at the center 
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of the MSB, and this shifts the lowest resonance, located 0.12 eV above the bottom of the well, to 
coincide roughly with the incident energy E * 0.005 eV. When this happens then resonant 
tunneling through the first barrier is highly probable, and there is a corresponding peak in the J-V 
characteristic. As the voltage is increased and the resonance is shifted below the incident energy 
level, resonant tunneling declines, giving rise to the region of negative differential resistance. The 
broad current peak at 0.86 V is likewise caused by the coincidence of the similarly broad 
resonance at 0.43 eV, with the incident energy; in other words, when the voltage dropped at the 
center of the MSB is 0.43 eV. 

It is worth pointing out, that these features of the J-V curve indicate unambiguously that 
the Stark shift of the MSB resonances is linear in the applied voltage or field strength. Also 
implied, is the fact that the shift, under localized fields, is simply caused by a shifted baseline or 
reference level. This is observed by others; for instance, by Schwartz (17), who reports a linear 
Stark shift in the resonances of AIGaAs MSB's, when calculating the resonances by means of a 
very similar transfer matrix technique. It appears likely, however, that if these calculations were to 
be repeated for MSB's subjected to uniform fields, the Stark shift of the resonances would not be 
linear, analogously to the behavior of bound states in single quantum wells. 

The J-V curves obtained by Esaki and Tsu are shown in Figure 47. Comparison with 
Figure 46 shows good agreement with the J-V curves calculated in the present study. This is a 
significant validation of the transfer matrix method used here. The curves of Esaki and Tsu show 
sharper peaks and deeper minima, however. This may be caused by two things. One is that they 
did not use a single incident energy, but integrated, according to equations 133a and b, the 
contribution of conduction electrons within about +/- kT of the Fermi energy. Another is that they 
neglected the voltage gradients inside the barriers and wells, taking account only of the drop 
between the centers of each. Although this is not a bad approximation for the barriers, which are 
only 2 nm wide, it may be less satisfactory for the 5 nm well regions. 

3.D.2. GaAs-AI(0.4)Ga(0.6)As MSB's and the Influence of the Effective 

Mass 

The resonances of the step barriers and the bound states of the single quantum wells 
discussed in the preceding sections, proved to be strongly influenced by the values of the 
effective mass used in calculating them. That the effective mass is similarly important to the 
outcome of tunneling calculations for MSB's, is shown by the T-E curve calculated for a two-barrier 
MSB, identical in geometry to the Esaki-Tsu GaAs-AI(0.5)Ga(0.5) As MSB, but with an effective 
mass m e ff which is constant, and equal to mo. These calculations are done in the program MSB. 
Figure 48 shows this T-E curve and the one belonging to the two-barrier GaAs-AI(0.5)Ga(0.5)As 
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MSB for which the proper effective masses are used. There is quite a significant difference 
between the two. Instead of two resonances, there are now five below the barrier height. The 
transmission coefficient is also lower by several orders of magnitude at most energies away from 
resonance. Analogously to quantum wells, the number of resonances that an MSB can support is 
proportional to the effective mass of the particle interacting with it. Similarly, the values of 
transmission minima through a single step barrier are a function of the diff erenc e between the 
effective masses inside and outside the barrier, and it appears that this is also the case for MSB's. 

The data in Figure 48 above, illustrate the consequences of using a grossly wrong 
effective mass in tunneling calculations for MSB’s. Suppose instead that the electronic effective 
mass of pure GaAs, m e ff - 0.0636, is used throughout, neglecting the change to m e ff = 0.0935 
inside the barriers. This is the approximation frequently made in the calculations of the bound 
states of single AIGaAs quantum wells, for instance in the work of Austin and Jaros and others (41 - 
43). It has already been shown that this strongly affects the location of the bound state in a single 
GaAs-AI(0.5)Ga(0.5) As quantum well. 

Figure 49a shows the T-E curves for two MSB's: one of them the two-barrier Esaki-Tsu 
MSB using the correct effective masses, and the other of the same geometry, but using m e ff - 
0.0636 everywhere inside and outside the barriers. The resonances in the latter case are more 
widely spaced and slightly broader than those of the realistically modelled MSB, with the lowest 
resonance shifted down to 0.10 eV, and the higher one to approximately 0.82 eV. The off- 
resonance transmission coefficient is also higher throughout, because the average effective 
mass is lower. 

Figure 49b shows the resulting J-V curves for the above two structures. The current 
peaks follow the same trends as the resonances of the T-E curves. When the effective mass is 
taken as 0.0636 everywhere, the low- voltage current peak is shifted to 0.20 V from 0.23 V. The 
tunneling current is approximately half an order of magnitude higher at all voltages, following the 
behavior of the transmission coefficient, and the current peaks are broadened in the same 
manner as the resonances in the T-E curve. 

The T-E and J-V curves of Figure 49 show that it is necessary to take careful account of 
the variations in effective mass in MSB's for accurate determination of the resonances. The MSB's 
used to illustrate this, however, are somewhat unrealistic, because the barrier heights are not 
correct for the assumed aluminum concentrations. These MSB’s served mainly as a basis for 
comparison with the results of Esaki and Tsu. Tunneling calculations for GaAs-AI(0.4)Ga(0.6)As 
MSB's are now discussed. These are layer compositions for which the maximum direct conduction 
band edge discontinuity, and hence, the largest possible barrier height, is 0.33 eV. Single step 
barriers of the same composition were discussed in section 3.B. above. The discussion to follow 
on such MSB’s is aimed at providing a basis for comparison with sawtooth barriers which have this 
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as the composition at the end of each layer. Results of tunneling calculations for these sawtooth 
superlattices are found in the next section. 

The T-E and J-V curves are calculated in the programs ALMSB and JVALMSB for GaAs- 
AI(0.4)Ga(0.6)As MSB's. These MSB's have two and five barriers, 2nm wide and 5 nm apart (i.e., 
with the same number, width and spacing as the Esaki-Tsu MSB's just dealt with), and barrier 
heights of 0.33 eV. These curves for the two-barrier MSB are presented in Figure 50a and b, 
respectively. Shown together with them, for the purpose of comparison, are the T-E and J-V 
curves for the two-barrier GaAs-AI(0.5)Ga(0.5)As MSB's. 

In Figure 50a, the upper T-E curve belongs to the GaAs-AI(0.4)Ga(0.6)As MSB. The 
resonances are shifted slightly to lower energies, and are broadened: the lowest resonance 
occurs at roughly 0.09 eV, and the upper one at 0.37 eV. The off-resonance transmission 
coefficient is also raised by at least an order of magnitude. In this case, these differences are the 
result not just of changing the value of the effective mass inside the barrier, (and hence also the 
average value for the MSB), but also of the fact that the barrier height is lower, at 0.33 eV instead 
of 0.5 eV. The lower average effective mass and the lower barrier height both contribute to 
greater probability of tunneling through this structure than in the Esaki- Tsu MSB. 

In Figure 50b, again the upper curve belongs to the GaAs-AI(0.4)Ga(0.6)As MSB. The 
features of this simple J-V curve follow from the T-E curve: there is a current peak at approximately 
0.18 V, or at twice the energy of the lowest resonance, and another very broad one around 0.75 
eV. The tunneling current is generally one half to one order of magnitude higher throughout than 
for the GaAs-AI(0.5)Ga(0.5)As MSB. 


3.E. Sawtooth Superlattices 

In this final section, the results of tunneling calculations for sawtooth barrier superiattices 
are presented. The transmission coefficient as a function of incident energy, and the tunneling 
current as a function of the applied voltage are computed, and the results are compared with 
those obtained for step barriers. 

The transmission coefficient is calculated by the matrix method outlined in section 
2.B.6.d., in the programs SAWTOOTH and ALSTAIR. These programs are found in Appendix A 
along with a description of their operation. Both programs accept the number of barriers, N, and 
assume right triangles, accepting their base and altitude. The applied electric field enters as a 
parameter. Data files of the transmission coefficient versus the incident energy are created and 
used to make T-E plots parametric in the applied field. The program SAWTOOTH assumes a 
constant effective mass equal to mo, the free electronic mass. The program ALSTAIR models 
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graded band-gap AI(x)Ga(1-x)As superlattices, and accepts the maximum and minimum aluminum 
concentrations within each barrier, assuming a linear concentration variation. The effective mass at 
each point in the barrier is calculated by means of equation 122. The modified connection rules 
taking account of the effective mass variation at each step are used in this program. 

The tunneling current is calculated as a function of the applied voltage in the programs 
JVSAWTOOTH and JVALSTAIR. These programs are derivatives of SAWTOOTH and ALSTAIR, 
respectively. Instead of the incident energy, the applied voltage is the independent variable. An 
incident energy of E = 0.005 eV is used just as for the MSB's of the preceding section. Data files 
of the tunneling current versus applied voltage are used to construct J-V plots. 

Below, results of tunneling calculations for single sawtooth barriers are first discussed and 
are compared with those already obtained for single square step barriers, in section 3.B. The 
effects of external electric fields and the effective mass are discussed. Then, tunneling 
calculations for multiple sawtooth barriers are presented, and compared to those for similar 
multiple step barriers, in particular the MSB's studied by Esaki and Tsu (11). 

3.E.I. Single Sawtooth Barriers 

3.E.1.a. Single Sawtooth Barrier In Zero External Field 

There is no analytical solution for the transmission coefficient through a sawtooth barrier, 
hence transfer matrices are especially well suited for this problem. It is interesting to compare the 
results of tunneling calculations for sawtooth barriers with those of step barriers, and to see what 
are the effects of the graded barrier height. Recall that in section 3.B., the transmission 
coefficient for 10 nm wide, 0.33 eV step barriers was calculated. The T-E curves obtained showed 
resonances above the barrier height: under zero field these occurred whenever the barrier width 
was an integral number of half wavelengths. The character of these curves was strongly 
influenced by the average effective mass as well as the difference in mass inside and outside the 
barrier. The behavior of waves encountering sawtooth barriers might be expected to be similar. 
The transmission coefficient for a lOnm wide, 0.33 eV sawtooth barrier with effective mass 
constant and equal to mo, is shown in Figure 51a, along with the T-E curve from Figure 28 for a 
step barrier of the same dimensions and effective mass. There is a striking absence of resonances 
in the T-E data of the sawtooth barrier. Tunneling is much more likely through the sawtooth barrier, 
for energies well below the peak barrier height and the maximum value of 1 is approached 
monotonically. By way of contrast, the transmission coefficient through the step barrier is 
extremely small until just above the barrier height, where it rises steeply to 1 , and then oscillates 
strongly. 
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The smooth shape of the transmission curve for the sawtooth can be thought of as 
resulting from the lack of a resonance condition such as exists in square step barriers. In fact a 
similar effect has already been encountered in step barriers under external, constant electric 
fields. There, the field effectively tapers the upper part of the barrier: the stronger the field, the 
broader the resonances and the lower the contrast between peaks and valleys. The approach to 
the maximum value of 1 is also more gradual as the field strength increases. The grade in a 
sawtooth barrier is analytically indistinguishable from the grade in a step barrier induced by an 
externally imposed electric field, so the similarity of their T-E data is not terribly surprising. 

Figure 51b shows the T-E data obtained for a single GaAs-AI(0.4)Ga(0.6)As sawtooth 
barrier, 10 nm wide and 0.33 eV high. The effective mass in this barrier is graded parabolically from 
0.0636 to 0.0871 at the peak. The curve for this barrier is superposed on that of a 1 0nm wide 
GaAs-AI(0.6)Ga(0.6)As step barrier. The same qualitative differences observed between the T-E 
data of the constant m e ff, step and sawtooth barriers in Figure 51a, are also present in Figure 51b. 
Considerably more tunneling takes place at all energies below the barrier height in the sawtooth, 
and there is a similar lack of resonances. 

In fact, as Figure 52 shows, tunneling through sawtooth barriers appears to be much less 
susceptible to the influence of the effective mass than is tunneling in step barriers. In Figure 52 
the T-E curves of the constant mass, m e ff = mo and the GaAs-AI(0.4)Ga(0.6)As sawtooth barriers 
are plotted together. The two curves are quite similar. Just as observed in step barriers, the lower 
average effective mass enhances tunneling at all energies. The transmission coefficient is 0.5 at 
0.25 eV in the GaAs-AI(0.6)Ga(0.6)As sawtooth, only reaching this value at 0.33 in the constant 
mass sawtooth. The similarity between the two curves is otherwise much stronger than between 
T-E curves belonging to comparable step barriers (Refer to Figure 30a). 

3.E.1.b. Single Sawtooth Barriers under Constant, Localized Electric 

Fields 

Field-assisted tunneling through sawtooth barriers is expected to differ somewhat from 
that through step barriers. If the field is constant, it results in an additional linear potential in the 
region of the sawtooth. This potential just adds to the already existing one, and as pointed out 
above, is not analytically different from it. This is different from imposing such a field on a step 
barrier, since that results in a situation that is analytically different from the zero-field case. For this 
reason it is expected that the effects of the field on tunneling through sawtooth barriers should be 
less radical than when dealing with step barriers. 

In this study, the applied field is taken as constant and opposing the gradient of the 
sawtooth itself. This is illustrated in Figure 2. As the field strength is increased to the value of the 


61 



sawtooth potential gradient, the barrier changes to a staircase potential. When this happens, 
tunneling no longer occurs, since there is no longer a barrier to tunnel through. The transmission 
coefficient must always be greater than zero for such a potential. Again, this is unlike field-assisted 
tunneling through step barriers, where no matter how great the field strength, a barrier to 
tunneling always exists. 

Field-assisted tunneling calculations for 10 nm, 0.33 eV sawtooth barriers are done for 
field strengths of 0.02 and 0.05 eV/nm. (These are same values used in the field-assisted 
tunneling calculations for the 10 nm, 0.33 eV step barriers.) Calculations are done for both m e ff = 
mo, and for the GaAs- AI(0.4)Ga(0.6)As sawtooth. The resulting T-E curves are shown for all field 
strengths in Figure 53. As the field strength is increased, the transmission coefficient also 
increases at all energies, but the approach to 1 is slowed. Note that the T-E curves for F = 0.05 
eV/nm are of a different character Than those obtained at 0.02 eV/nm: T is greater than zero for all 
energies in the former case. Figure 54 shows the deformation of the sawtooth barrier under the 
applied field. At F = 0.05 eV/nm, the structure is a staircase, not a sawtooth, this transition occurs 
at a field strength of 0.033 eV/nm. At field strengths less than this, tunneling can still occur, and, 
for the constant mass sawtooth, the T-E curves have an inflection point, at energies equal to the 
peak height of the sawtooth, where the value of the transmission coefficient is 6 . 5 . For the GaAs- 
A((0.4)Ga(0.6)As sawtooth, these inflection points are shifted to lower energies. 

An interesting feature of the curves in Figure 53, is that as the field strength increases, 
the difference between the T-E curves of the GaAs- AI(0.4)Ga(0.8)As and the constant mass 
sawtooth barriers becomes less significant. The curves obtained at F « 0.05 eV/nm in fact are 
practically identical. The reason for this is not apparent. This trend was not evident in the step 
barriers studied here, although it might have been, in the limit of very high fields, which was not 
adressed in this study. 

Overall, it appears that single sawtooth barriers differ from step barriers in a few respects. 
One is that, because of the sawtooth's graded thickness, transmission resonances are almost 
totally supressed. This is a consequence of the lack of a resonance condition such as exists in 
step barriers. Another is that the interaction of the incident wave with the sawtooth barrier 
changes character at field strengths sufficient to deform the barrier into a staircase structure. 
Finally, for reasons that are not clear, the influence of the effective mass is very much less 
pronounced, and tends to become still less important as the field strength is raised. In the next 
section, the same issues will be adressed for multiple sawtooth barriers. 


3.E.2. Multiple Sawtooth Barriers 
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Recalling that resonances below the barrier height and equal to one were possible in 
superlattices of step barriers, one expects to observe a similar phenomenon in sawtooth 
superlattices. The resonance peaks, however, are anticipated to be broader than those in MSB's, 
because instead of square wells there are now triangular ones. Triangular barrier shapes result in 
suppressed and broadened transmission resonances as was observed in single sawtooth barriers 
and step barriers under external fields. Thus any resonance peaks found in sawtooth 
superlattices are expected to be broader than those in similar MSB's of the same composition. 

In Figure 55, the transmission coefficient as a function of energy is plotted for a two-barrier 
sawtooth superlattice, with bases equal to 4.5 nm, peak heights of 0.5 eV, and m e ff taken as 
constant and equal to mg. Compare this with the same data for a two-barrier MSB, with 2nm wide 
bases, 5nm apart, barrier heights of 0.5 eV, and m e ff equal to mo, found in Figure 48. There are 
indeed resonances in the transmission coefficient of the sawtooth superlattice, but they are fewer 
in number than those of the MSB. Instead of five, there are only two occurring below the barrier 
height of 0.5 eV. The upper resonances are considerably broader than those of the MSB. Finally, 
the value of the transmission coefficient is higher, away from the resonances, than in the MSB, 
except at the lowest energies, where the thickness of the sawtooth and the step barriers is the 
same. These discrepancies between the two curves in Figure 55 are explained by the difference 
in shape of the sawtooth and step barriers. 

In Figure 56, the T-E data for two-barrier sawtooth superlattices with differing effective 
masses are plotted. The data shown are for the same geometry as above, with Figure 56a 
showing the transmission coefficient for m e ff * mg, and Figure 56b showing that of a superlattice 
with the composition in each barrier graded from pure GaAs to AI(0.5)Ga(0.5)As. The effective 
mass is thus graded from 0.0636 to 0.0935 in each sawtooth. Comparing Figures 56a and b 
shows that the effective mass variation eliminates two of the resonances occurring below the 
barrier height. This is most likely a result of simultaneous shifting and broadening of the 
resonances. A single, very broad resonance at approximately 0.31 eV persists, however. 

It appears that in AI(x)Ga(1-x)As sawtooth superlattices, the effective mass variations have 
a rather more pronounced effect on tunneling than in single sawtooth barriers. As Figure 53 
illustrated, the T-E curves obtained for single sawteeth were quite insensitive to the values used 
for the mass. This is evidently not the case for tunneling through multiple sawtooth barriers. 

3.E.2.a. GaAs-AI(0.5)Ga(0.5)As Sawtooth Superlattices: Comparison with 

Similar MSB's (Esakl-Tsu) 

Comparison is made in this section between G?As-AI(0.5)Ga(0.5)As step and sawtooth 
superlattices. Since the tunneling properties of such MSB's are well understood, and calculations 
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have been done for them here and elsewhere, they provide a convenient basis for comparison 
with sawtooth superlattices. 

The transmission coefficient is calculated for two- and five-barrier sawtooth superlattices, 
with bases of 4.5 nm, heights of 0.5 eV, and with composition graded from GaAs to 
AI(0.5)Ga(0.5)As. These dimensions result in the same number of periods and overall length for 
sawtooth superlattices, as for MSB's with 2nm wide barriers 5nm apart. The transmission 
resonances and consequent tunneling current-voltage characteristics should thus be somewhat 
similar. 

The T-E curves for the two-barrier step and sawtooth superlattices are shown in Figure 
57. Instead of the two resonances supported by the MSB, the sawtooth superlattice has only a 
single very broad one around 0.31 eV. The effective mass gradient, as well as the barrier gradient, 
likely contribute to the breadth of this resonance. 

The T-E curves for the five-barrier GaAs-AI(0.5)Ga(0.5)As sawtooth superlattice and MSB 
are found in Figure 58. The MSB supports two sets of resonances below the barrier height, each 
resonance split by coupling into four peaks. These are centered around 0.12 and 0.43 eV. The 
sawtooth superlattice, however, has only a single set of resonances. These are likewise four in 
number and are centered around the two-barrier resonance at 0.31 eV. These four resonances 
are much broader than those of the MSB. The transmission coefficient is higher for the sawtooth 
barrier down to an energy of approximately 0.22 ev, below which energy the sawteeth are thicker 
than the step barriers. 

The J-V curves for the two-barrier structures are found in Figure 59. These curves are 
both calculated assuming an electric field localized to the superlattice. Note that there is a 
qualitative similarity between the two curves: a low-voltage current peak, followed by an area of 
negative differential resistance (NDR), then another broader peak at higher voltage, also followed 
by a regime of NDR. The current peaks of the sawtooth J-V curve are naturally broader than those 
of the MSB, since the resonances are very broad. The appearance of regions of NDR in sawtooth 
superlattices means that they may find novel applications in electronics, just as MSB's have (6). 

It was shown earlier that the features of the J-V characteristic of this MSB were explained 
by resonant tunneling between the incident energy level (0.005 eV), and the resonances shifted 
downward by the applied field. The fact that the voltage peaks occur at twice the resonance 
energy levels, implied a linear shift of the resonances under the field. 

Comparison of the J-V characteristic of the sawtooth superlattice with its T-E curve, 
however, indicates a different behavior under the field. There is a current peak at about 0.47 volts, 
and another at roughly 1 .3 volts. Neither of these would result from a resonance at 0.35 eV 
shifted linearly downward with the applied field. It appears that for sawtooth superlattices, as the 
applied field strength is raised, new resonances form in the triangular wells, which may, at certain 
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values of the field, coincide with the incident energy of 0.005. In other words, the Stark effect in 
sawtooth superlattices is more complicated than in MSB's. This is a consequence of the way the 
sawtooth deforms under external fields. As the field strength is raised, the shape of the well 
between the barriers changes; In particular, the slope of the leftmost edges decreases. This is 
unlike the behavior of MSB's under uniform fields, as Rgure 3 shows. 

Figure 60 shows the shape of this sawtooth superiattice under applied voltages of zero, 
0.47, and 1 .3 V. At 0.47 V, there is still a small barrier to tunneling, but at 1 .3 V the structure is a 
staircase and no tunneling is occurring. At 0.47 V, a shallow state is still probably trapped between 
the barriers, from which resonant tunneling can take place. 

Figure 61 shows the T-E curve for this sawtooth superiattice under an an applied field 
strength of 0.052 eV: this corresponds to an applied voltage of 0.47 V. A slope change is evident 
at approximately 0.03 eV; another at around 0.35 eV. A shallow bound state may be hinted at by 
the knee in the curve at 0.03 eV. If such a state does exist, however, it must necessarily be 
difficult to resolve in the T-E data, since its energy level must be close to zero ( i.e.to the incident 
energy level of 0.005 eV). 

The peak at 1 .3 V is a bit more puzzling, since no trapped states exist at such a high field. 
As earlier observed in single step barriers, however, there can still be pronounced resonances at 
energies well above the barrier height. This may explain the broad peak and region of NDR at 
around 1 .3 V. 

3.E.2.b. GaAs-AI(0.4)Ga(0.6)As Sawtooth Superlattices 

The preceding section explored the similarities between sawtooth superlattices and 
multiple step barriers. It focused on GaAs-AI(0.5)Ga(0.5)As superlattices mainly because of the 
abundance of existing data on MSB's of that composition with which comparisons could be made. 
AI(x)Ga(1-x)As is an indirect gap material for aluminum concentrations greater than approximately 
0.4 mole fraction.however.and only tunneling between direct conduction bands is considered in 
this work. Hence this final section will focus on GaAs-AI(0.4)Ga(0.6)As sawtooth superlattices, as 
this is the composition giving the greatest direct conduction band discontinuity of about 0.33 eV. 

The transmission coefficient of a two-barrier sawtooth superiattice, with 4.5 nm bases and 
peak heights of 0.33 eV, is shown as a function of energy in Figure 62. Superposed is the T-E 
cun/e of a similar GaAs-AI(0.5)Ga(0.5)As MSB. As might be expected, the curves are qualitatively 
very similar. The resonance in the AI(0.4)Ga(0.6)As data is shifted to approximately 0.23 eV, 
however, compared to 0.31 eV for the GaAs-AI(0.5)Ga(0.5)As MSB. Off resonance, the 
transmission coefficient is also higher throughout the range. These changes are consequences 
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of the lower average effective mass throughout each barrier, as well as of the lower barrier height, 
in GaAs-AI(0.4)Ga(0.6)As. 

The tunneling current as a function of the applied voltage in the two- barrier superlattices 
is shown in Figure 64. Again, a strong qualitative resemblance between the two curves is 
observed. The tunneling current is greater in the AI(0.4)Ga(0.6)As superlattice by about two 
orders of magnitude, and this naturally follows from the higher transmission coeffident observed 
in Figure 62. A low-voltage current peak at roughly 0.35 V, followed by a region of negative 
differential resistance is evident in the AI(0.4)Ga(0.6)As curve. As was the case for the 
AI(0.5)Ga(0.5)As sawtooth superlattice, the location of this peak implies a more complicated 
change In the resonances under applied fields than in MSB's. Certainly the shift is not a simple 
linear function of the field strength. 
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4. CONCLUSION 


A transfer matrix method has been applied to quantum mechanical tunneling calculations 
in GaAs-AI(x)Ga(1 -x)As semiconductor heterostructures. The focus of this work was a novel type 
of heterostructure in which the band gap is graded linearly, giving rise to a sawtooth superlattice. 
Tunneling calculations for such a structure are presented here for the first time. In addition, 
tunneling calculations were performed for conventional step barrier heterostructures. The Stark 
shift of the bound and quasi-bound states of single finite quantum wells was calculated as well. 
Some conclusions are: 

1 . The Stark shift of bound states in single quantum wells can be obtained by the transfer 
matrix method used here. The extent of the applied field, however, determines the functional 
dependence of the Stark shift. If the field is restricted to an area near the well, the Stark shift is 
linear in the field. If the field is infinite in extent, the shift is quadratic. The extent of the field must 
be considered in tunneling calculations for sawtooth and step barrier superlattices as well. 

2. The tunneling properties of step and sawtooth superlattices show some strong 
qualitative similarities. Both structures can exhibit resonant tunneling, as evidenced by correlation 
of transmission resonances with peaks in the current-voltage curves. However, the shift of the 
resonances is linear in the field in step barrier superlattices, while in sawtooth superlattices the 
shift is not a simple function of the field. This is because the two kinds of structure deform 
differently under uniform fields. The sawtooth deforms into a staircase under a high enough field 
strength, and tunneling no longer occurs since the barriers are eradicated. Step barriers always 
present some barrier to tunneling no matter how strong the field. 

3. The effective mass variations encountered in semiconductor heterostructures should 
not be neglected. The conventional wave function boundary conditions at interfaces must be 
modified to conserve probability current density when the mass is discontinuous. The range of 
effective mass in the GaAs-AI(x)Ga(1-x)As heterostructures studied here was found to have a 
significant effect on the outcome of calculations. 

a. The number and level of the bound states in single quantum wells is strongly affected 
by average effective mass, and by the discontinuity at well edges. Neglecting the average 
effective mass results in a gross error in the number of levels, while ignoring discontinuities at the 
well edges results in a significant change in their energy. 

b. The use of the correct effective masses in single and multiple step and sawtooth 
barriers is also important ,as the number and energy of transmission resonances is a function of 
the average effective mass and the discontinuities at interfaces In sawtooth superlattices the 
parabolic effective mass gradient within each layer should be taken account of. 
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1. Quantum Confinement and Resultant Quantization of Energy Eigenvalues (1) 
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3. Step- and Sawtooth-Barrier Superlattices under Reverse Bias 

a. Step-barrier superlattice under zero external field. Virtual states are shown; 

b. Step-barriers under strong reverse bias. Virtual states can still exist; 

c. Sawtooth-barrier superlattice under zero field. Virtual states are formed as in the case of step- 

barrier superlattices; 

d. Sawtooth-barrier superlattice under moderate reverse bias. The applied voltage is less than the 

sum of the barrier heights, and virtual, or quasi- bound, states can still be formed. 

e. Sawtooth-barrier superlattice under strong reverse bias. Structure is actually a staircase and 

quasi-bound states no longer are possible, although resonances still occur. 
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4. Transfer Matrix at a Potential Step 
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a. Transmission coefficient for uniform effective masses; 

b. Transmission coefficient for different effective masses inside and outside the barrier. 


(-i/kj 
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V(z) 
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quasi-bound, states, while MQWs can have true bound states. Both types of heterostructure can 
support resonances for E greater than the barrier or well height. 



Criterion for bound states 
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11. Bound States and Resonances of a Single Quantum Well. 

a. E > 0: Transmission coefficients for uniform and varying effective masses. Resonances in T ; 

b. E < 0: Criterion for bound states. Eigenvalue condition on Mu gives bound states. 
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The region inside the well is divided into intervals where plane-wave transfer matrices are 
calculated. In the regions outside the well, the solutions of the Schroedinger equation are the Airy 
functions. 
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18. Deformation of a Single Quantum Well In a Localized Electric Field. 

The applied field shifts the bound and quasi-bound levels to lower energy. 

a. Low field strength, with true bound states at Eo and Ei; 

b. Moderate field strength, with one true bound state at Eo - A, and a quasi-bound state at Ei - A; 

c. High field strength, only quasi-bound states remain. 
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External, Localized Electric Field F 






22. Energy Band Structure of Pure GaAs (<100> and <111> directions) (33). 
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(a) (b) 

23. Complex Energy Band Structures of pure GaAs and AlAs (34) 

i. (110) interface: a. GaAs b. AlAs; 

ii. (100) interface: a. GaAs b. AlAs. 
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24. Energy Gap In AI(x)Ga(1-x)As as a Function of AlAs Mole Fraction x (19). 

The x-dependence of the direct conduction band Ti c is shown by the solid line; that of the 
indirect gap X-| C by the dashed line. The direct and indirect minima are equal at x = 0.37. 
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25. The Fraction of Conduction Electrons In r 1c of AI(x)Ga(l-x)As as a Function 
of AlAs Mole Fraction x.( 18 ) 

Data are taken at 300 K. Dotted line is for degenerate case with N = 4xi0 17 cm' 3 ; solid line is for 
nondegenerate case with N = 4x1 o 1 6 cm' 3 . 
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ACAs mole fraction 


26. Band-Edge Alignments at GaAs-AI(x)Ga(1-x)As Heterojunctions (R. Miller, 
ATT-Bell Laboratories). 

a. AE C , conduction band misalignment; 

b. AE V , valence band misalignment. 
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27. Esakl*Tsu Multiple Step Barrier Geometry (11) 

a. Zero applied electric field; 

b. Applied field strength F = F ap over the length I of the MSB. 
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28. T-E Data for Single Step Barrier. 

Barrier is 10 nm wide, 0.33 eV high. Effective mass is uniform and equal to free electronic mass 
mo- 
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29. T-E Data for Single GaAs-AI(0.4)Ga(0.6)As Step Barrier. 

Barrier is 10 nm wide and 0.33eV high. Effective mass is mj n = 0.0871 inside the barrier, rriout = 

0.0636 outside. 
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(b) E (eV) 


30. T-E Data for Single Step Barriers lOnm Wide and 0.33 eV High. 

a. Superposed T-E curves for GaAs-AI(0.4)Ga(0.6)As barrier, and barrier with uniform effective 
mass m e ff = mo everywhere. Curve 1: GaAs-AI(0.4)Ga(0.6)As, mj n = 0.0871, m^t = 0.0636; 
Curve 0: m e ff = mo; 

b. Superposed T-E curves for GaAs-AI(0.4)Ga(0.6)As barrier, and barrier with uniform effective 
mass m 0 ff = 0.0636 everywhere. Curve 0: GaAs-AI(0.4)Ga(0.6)As; Curve 1 : m e ff = 0.0636. 
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31. Effects of Applied Electric Field on Transmission Coefficient of Single Step 
Barrier. 

Barrier is 10 nm wide, 0.33 eV high. Effective mass = mo everywhere. 

a. Fap = 0; 

b. F ap = 2 xlO' 2 eV/nm; 

c. Fap = 5 xl O' 2 eV/nm. 
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32. Effects of Applied Electric Field on Transmission Coefficient of Single 
GaAs-AI(0.4)Ga(0.6)As Step Barrier. 

Barrier is 10 nm wide, 0.33 eV high. Inside the barrier m^f = 0.0871, outside the barrier m e ff = 
0.0636. Curve 0: F a p = 0; Curve 1 : F a p = 2 xIO" 2 eV/nm; Curve 2: = 5 xl O' 2 eV/nm. 
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33. Applied Field and Effective Mass Effects on Transmission Coefficient of 
Single Step Barriers. 

Barriers are lOnm wide, 0.33 eV high. The two curves in each figure are for a barrier with uniform 
effective mass m e ff = mo. and a GaAs-AI(0.4)Ga(0.6)As barrier: Curve 0: m e ff = mo; Curve 1 : mj n = 
0.0871, rriout = 0.0636. 

a. Applied field F ap = 2x1 0‘ 2 eV/nm ; 

b. Applied field F ap = 5 xlO' 2 eV/nm . 
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34. Total Transfer Matrix Element M^ for Single GaAs-AI(0.5)Ga(0.5)As 
Quantum Well: Bound State. 

Well is 3 nm wide, 0.4 eV deep. This calculation took account of the effective mass variation at the 
well edges: the effective masses inside and outside the well are mj n = 0.094, m out = 0.0636. The 
minimum of Mi i gives a bound state energy of -0.197 eV below the top of the well. 
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35. Bound State Energy as a Function of the Effective Mass Inside Well. 

Well is the same as that of Figure 34. The effective mass, however, is taken in this calculation as 
0.0636 inside and outside the well. The new bound state energy is -0.217 eV below the top of 
the well, instead of -0.197 eV. 

Curve 0: Effective mass = 0.0636 everywhere. E 0 = -0.217 eV. 

Curve 1: Effective mass = 0.094 inside the well, 0.0636 outside the well.Eo = -0.197 eV. 
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36. Bound States In Single Quantum Well. 

Well is 3nm wide, 0.4 eV deep. Effective mass is uniform and equal to mo. The minima in Mi i 
correspond to bound states. Using the free-electronic mass results in four bound states, in 
agreement with Equation 139. 
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37. Total Transfer Matrix Element Mu as a Function of Localized Field Strength. 
Data are for single GaAs-AI(0.5)Ga(0.5)As quantum well, 3 nm wide and 0.4 eV deep. Effective 
mass taken as 0.0636 everywhere. Each curve corresponds to a different field strength. Minima 
shift to energies deeper in the well as the field strength is increased. 

Curve 0: Fgp = 0; Curve 1 : F ap = 1 xl O' 2 eV/nm; Curve 2: Fgp = 2 xIO' 2 eV/nm; 

Curve 3: F ap = 3 xIO' 2 eV/nm; Curve 4: Fgp = 4 xIO' 2 eV/nm; 

Curve 5: F ap = 5 xIO' 2 eV/nm. 
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38. Stark Shift of Bound State Energy Level under Localized Field. 

Data are for single GaAs-AI(0.5)Ga(0.5)As quantum well, 3 nm wide and 0.4 eV deep. Effective 
mass taken as 0.0636 everywhere. In this figure the Stark shift (referred to the zero-field bound 
state energy) is plotted as a function of the applied field strength. The Stark shift is linear in the 
field strength. 
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39. The Ratio (F/E) as a Function of Energy and Applied Field Strength. 

Data are for single GaAs-AI(0.5)Ga(0.5)As quantum well, 3 nm wide and 0.4 eV deep. Effective 
mass taken as 0.0636 everywhere; field is uniform and unrestricted. Each curve is for a different 
field strength; the maxima in the data correspond to quasi-bound states of the quantum well. 
Maxima shift to energies deeper in the well as the field strength increases. 

Cun/e 0; F ap = 1.0 xl O' 2 eV/nm; Cun/el: F ap =1.5 xl O' 2 eV/nm; 

Curve 2; F ap = 2.0 xIO' 2 eV/nm; Curve 3: F ap = 2.5 xl O' 2 eV/nm; and 
Curve 4; F ap = 3.0 xl O’ 2 eV/nm . 
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40. Stark Shift of Bound State Energy Level Under Uniform, Unrestricted Field. 
From Figure 39, the energy levels of the maxima in (F/E) are plotted as function of the applied field 
strength. The bound state energies are quadratically dependent on the applied field strength. 
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41. Stark Shift of Bound State Energy Levels Under Uniform, Unrestricted Field. 

From Figure 40, the Stark shift relative to the zero-field level is plotted as a function of the applied 
field strength. Also shown are the data of Austin and Jaros (41) for an identical quantum well. 
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42. Stark Shifts under Localized and Unrestricted Uniform Fields Compared. 
From Figures 38 and 41 , the Stark shifts calculated for a 3nm, 0.4 eV AI(0.5)Ga(0.5)As quantum 
well are plotted together. 
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43. T-E Curves for Two-, Three-, and Five-Barrier GaAs-AI(0.5)Ga(0.5)As 
Multiple Step Barriers. 

These are the results of transfer matrix calculations for step barrier superlattices identical to those 
of Esaki and Tsu (11). Barriers are 2 nm wide, 5 nm apart, and 0.5 eV high. Effective mass in the 
barriers is mj n = 0.094, between the barriers mo Ut = 0.0636. Modified connection rules are used. 
Figures a, b, and c are for two, three and five barriers, respectively. 
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44. T-E curves for Two* and Five- Barrier Esaki-Tsu Type MSB's 

The T-E curves from Figure 43 for two and five barriers superposed. Note the splitting of single 
resonances into four, caused by coupling between wells. 
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45. T-E Curves for Two-, Three-, and Five-Barrier MSB’s Calculated By Esaki 
and Tsu (11). 

Note the close agreement between these curves and those of Figure 43. 
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46. J-V Curves for Two- and Three-Barrier GaAs-AI(0.5)Ga(0.5)As MSB's 
These curves are calculated for the same MSB's whose T-E curves appear in Figure 43. 
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47. J-V Curves for Two- and Three-Barrier MSB's Calculated by Esaki and Tsu 
( 11 ). 

Note the agreement between these curves and those of Figure 46. 
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48. T-E Curves of Two-Barrier MSB's: Influence of Effective Mass. 

Both curves are for two-barrier MSB's with barriers 2 nm wide, 5 nm apart, and 0.5 eV high. Curve 0 
is for GaAs-AI(0.5)Ga(0.5)As MSB with effective mass variations taken into account, and Curve 1 is 
for MSB with effective mass uniform and equal to mo. 
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49. T-E and J-V Curves of Two-Barrier MSB's: Effect of Neglecting the Effective 
Mass Step at Heterojunctions. 

Data are for GaAs-AI(0.5)Ga(0.5)As Esaki-Tsu-type MSB's. Figures a and b are T-E and J-V curves 
respectively. In Curve 0, the effective mass is taken as uniform and equal to 0.0636. In Curve 1, 
the effective mass is 0.094 inside the barriers, and 0.0636 in the wells. 
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50. T-E and J-V Curves for GaAs-AI(0.4)Ga(0.6)As and GaAs-AI(0.5)Ga(0.5)As 
Two-Barrier MSB's: Effects of Composition 

Data are for MSB's with barriers 2 nm wide, 5 nm apart. Figures a and b are T-E and J-V curves 
respectively. Curve 0 is for GaAs-AI(0.4)Ga(0.6)As. Curve 1 is GaAs-AI(0.5)Ga(0.5)As. 
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51. T-E Curves for Sawtooth and Step Barriers Compared. 

Sawtooth barrier is 10 nm wide at the base, step barrier is 10 nm wide. Both are 0.33 eV high. 
Figure a. Eftective mass is uniform and equal to mo. Curve 1 is for the sawtooth, Curve 2 is for the 
step barrier. 

Figure b. Barriers are made of GaAs-AI(O.4)Ga(0.6)As. Effective mass is 0.0636 in GaAs, 0.0871 in 
AI(0.4)Ga(0.6)As. Curve 0 is for step barrier, Curve 1 is for the sawtooth. 
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52. T-E Curves for GaAs-AI(0.4)Ga(0.6)As Sawtooth Barriers: Influence of the 
Effective Mass. 

Both barriers are 10 nm wide at the base and are 0.33 eV high. Curve 0 is for GaAs- 
AI(0.4)Ga(0.6)As with effective mass variations taken account of, and Curve 1 is for effective mass 
uniform and equal to mg. Note the similarity between the two curves, as opposed to the same data 
for step barriers. 
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53. T-E Curves for Sawtooth Single Step Barriers: Influence of Effective Mass 
and Applied Electric Field Strength. 

All data is tor barriers 10 nm wide at the base and 0.33 eV high. Curve 0: F ap = 0, m e ff = mo- 
Curve 2: F ap = 0, GaAs-AI(0.4)Ga(0.6)As Curve 3: F ap = 0.02 eV/nm, m e ff = mo. Curve 4: F ap = 
0.02 eV/nm, GaAs-AI(0.4)Ga(0.6)As Curve 5: F ap = 0.05 eV/nm, m e ff = mo- Curve 6: F ap = 0.05 
eV/nm, GaAs-AI(0.4)Ga(0.6)As 
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54. Deformation of Single Sawtooth Barrier under Applied Electric Field. 

The shape of the barrier under F ap = 0, 0.02, and 0.05 eV/nm is shown in Figures a, b, and c, 
respectively. At Fap = 0.05 eV/nm the barrier is actually a staircase. 
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57 . T-E Curves for Sawtooth and Step Two-Barrier GaAs-AI(0.5)Ga(0.5)As 
Superiattices. 

Curve 0: MSB with barriers 2 wide, 5 nm apart, and 0.5 eV high 
Curve 1 : Sawtooth superlattice with 4.5 nm wide bases. 
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58. T-E Curves for Sawtooth and Step Five-Barrier GaAs-AI(0.5)Ga(0.5)As 
Superlattices. 

Curve 0: Sawtooth superlattice with 4.5 nm wide bases. 

Curve 1 : MSB with barriers 2 wide, 5 nm apart, and 0.5 eV high 
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59. J-V Curves for Two-Barrier GaAs-AI(0.5)Ga(0.5)As Sawtooth and Step 
Superlattices. 

Curve 0: Sawtooth superlattice with 4.5 nm bases, 0.5 eV high 
Curve 1 : MSB with barriers 2 nm wide, 5 nm apart, 0.5 eV high 
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60. Deformation of Two-Barrier Sawtooth Superlattice under Applied Electric 
Field. 

In Figures a, b, and c the applied field strength Fap is 0, 0.02 and 0.05 eV/nm, respectively. At F ap 
- 0.05 eV/nm the structure is a staircase. 
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61. T-E Curve for Two-Barrier GaAs-AI(0.5)Ga(0.5)As Superlattice 

Bases are 4.5 nm wide, heights are 0.5 eV. The applied field strength is 0.052 eV/nm (applied 

voltage = 0.47 V). 
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62. T-E Curves for Two-Barrier Sawtooth Superlattices: Effects of Composition. 
Data for GaAs-AI(0.4)Ga(0.6)As and GaAs-A!(0.5)Ga(0.5)As superlattices are shown in Curves 0 
and 1 respectively. Bases are 4.5 nm wide. 
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63. J- 


and 1, 
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APPLIED VOLTAGE (V) 

wo-Barrler Sawtooth Superlattices: Effects of Composition. 

(0.6)As and GaAs-Al(0.5)Ga(0.5)As superlattices are shown in Curves 0 
s are 4.5 nm wide. 
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APPENDIX 


The tunneling calculations of the present study use the programs in this appendix. All of 
the programs are based on the transfer matrix method described in Section 2.B.; their structure 
actually differs very little from program to program. These programs are written in VAX FORTRAN, 
and are compiled using the G_FLOAT option, which allows computations involving numbers 
roughly as large as 10®®®. This proves necessary in programs calculating resonance spectra and 
tunneling current voitage curves for superlattices longer than 10 nm or so. 

The transmission resonance programs are: ZIGGURAT; ZIGGEFFMASS; MSB; ALMSB; 
SAWTOOTH, and ALSTAIR. To calculate transmission resonance spectra, the incident energy is 
varied in small steps in a do-loop. The shape of the potential is determined from input parameters 
of barrier composition, number, height and width, and from the applied electric field strength. The 
potential is divided into many small steps of constant potential and effective mass. The transfer 
matrix, denoted by T in the code, is initialized as the unit matrix. In the main program, the position 
zj, wave numbers kj and kj+i , and effective masses mj and mj+i , are calculated at each step. 
These values are passed to a subroutine, EFFMASSTEP (or STPTRN if the effective mass is 
taken as uniform and equal to mo), where the transfer matrix at the step is computed, multiplied by 
the product of the transfer matrices at all preceding steps, and the accumulated T is returned to 
the main program. After this has been done at each step In the superlattice, the total accumulated 
transfer matrix element T-n is used to calculate the transmission coefficient (equation 18). This 
value is deposited with the corresponding incident energy E in an output file from which T-E plots 
are made. 

Tunneling current-voltage programs are: JVMSB, JVALMSB, JVSAWTOOTH, and 
JVALSTAIR. To calculate J-V curves, a single incident energy is chosen, 0.005 eV in this study. 
The applied voltage over the structure is varied within a do-loop in steps of 0.01 V, from 2 V to 
zero. Inside the do-loop, the operation of the program is very similar to that of the T-E programs: 
the potential shape is calculated from input parameters and divided into steps, the transfer matrix 
calculated at each step, inside the subroutines EFFMASSTEP or STPTRN, and the accumulated 
transfer matrix returned to the main program. After the last step, the tunneling current is calculated 
from Tii according to equation 20, and deposited with the corresponding value of V in an output 
file for plotting. 

The program WELL is used to calculate the bound states in a single quantum well where 
the effective mass is uniform and equal to mo. The programs for calculating the bound states in 
single GaAs-Al x Ga(i- x )As quantum wells are ALGAASWELL and STARKMASS. WELL and 
ALGAASWELL are used to calculate the bound states under zero applied field, while 
STARKMASS is used to find the Stark shift of those states under non-zero, localized electric 
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fields. These programs are almost identical to those used to calculate T-E curves for step barriers. 
The differences are that the energies are negative, i.e. Vo<E<0; then |togioCril)l (where Tn is 
from the total transfer matrix accumulated inside and outside the well) is simply deposited with it's 
corresponding energy in an output file. 

To calculate the Stark shift in a GaAs-Al x Ga(i- X )As quantum well under an extended 
electric field, the program AIRYWELL is used. In this program the Airy function solutions of the 
Schroedinger equation inside and outside the well are calculated from the input values of the well 
width and depth and applied field strength. The values of x, the argument of the Airy functions, 
defined by equation 72, are calculated at the well edges. The Airy functions are obtained from the 
IMSL, Incorporated, special function library SFUN/LIBRARY. Matching wavefunctions and 
derivatives at the well edges gives the wavefunction coefficients to the left of (C and D) and inside 
(P and Q) the well. From these and from R, S, T and U, the ratios of Airy functions evaluated at the 
well edges, the components F and E are calculated from the relation: 

F/E - (C/D * (S-P) + T - Q) / (U * (C/D * P + Q) - R * (C/D * S + T)) 

This is equivalent to equation 86 in the text. This calculation is done inside a do-loop which varies 
the energy in steps equal to 0.001 times the zero-field well depth plus the applied voltage across 
the well. Then |logio( F/E )l and the corresponding value of energy are stored in an array. In this 
program an NCAR graphics subroutine is used to obtain plots from the data array. 
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SUBROUTINE STPTRN { E, VL, VR, FK, Z, Tl) 
c 

c This subroutine computes the accumulated transfer matrix TR 

c through a potential step, where the wavenumber changes from AKL 

c on the left, to AKR on the right, at position Z. 

c 

IMPLICIT REAL* 8 (A-H, 0-Z) 

DIMENSION Tl (4) , T2(4), TR(4), EX(4) 

COMPLEX* 16 AKL, AKR, Tl, T2, TR, EX 
PSILON = 10E-7 
DIFFL - (E - VL) 

DIFFR - (E - VR) 

IF (ABS (DIFFL) .LT. PSILON) DIFFL * SIGN (PSILON, DIFFL) 

IF (ABS (DIFFR) . LT . PSILON) DIFFR - SIGN (PSILON, DIFFR) 

AKL - FK*DIFFL 
AKL = SQRT(AKL) 

AKR = FK*DIFFR 
AKR = SQRT(AKR) 

EX (1) = (0,1.0) * (AKR - AKL) * Z 

EX (2) - (0, -1.0) * (AKR - AKL) * Z 

EX (3) - (0, -1.0) * (AKR + AKL) * Z 

EX (4) - (0,1.0) * (AKR + AKL) * Z 

T2 (1 ) = 0.5 * (1 + AKR/AKL) * EXP(EX(1)) 

T2(2) - 0.5 * (1 + AKR/AKL) * EXP (EX (2)) 

T2 (3) = 0.5 * (1 - AKR/AKL) * EXP (EX (3)) 

T2 (4) = 0.5 * (1 - AKR/AKL) * EXP(EX(4)) 

TR(1) = Tl (1) *T2 (1) + Tl (3) *T2 (4) 

TR (2 ) - Tl (4 ) *T2 (3) + T1(2)*T2(2) 

TR (3) » Tl (1) *T2 (3) + Tl (3) *T2 (2) 

TR ( 4 ) - Tl (4) *T2 (1) + Tl(2)*T2(4) 

Tl (1) « TR(1) 

Tl (2) - TR (2 ) 

Tl (3) - TR(3) 

Tl (4) = TR(4) 

400 RETURN 

END 
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c 

c 

c 

c 

c 


400 


SUBROUTINE EFFMASSTEP ( E, VL, VR, ML, MR, FK, Z, Tl) 

This subroutine computes the accumulated transfer matrix TR 
through a potential step, where the wavenumber changes from AKL 
on the left, to AKR on the right, at position Z. 


IMPLICIT REAL* 8 (A-H, O-Z) 

DIMENSION Tl (4) , T2(4), TR(4), EX(4) 

COMPLEX* 16 AKL, AKR, Tl, T2, TR, EX 
REAL* 8 ML, MR 
PSILON = 10E-7 
DIFFL - (E - VL) 

DIFFR - (E - VR) 

IF (ABS (DIFFL) . LT . PSILON) DIFFL = SIGN (PSILON, DIFFL) 
IF (ABS (DIFFR) . LT . PSILON) DIFFR “ SIGN (PSILON, DIFFR) 
AKL - FK * ML * DIFFL 
AKL - SQRT(AKL) 

AKR = FK * MR * DIFFR 
AKR = SQRT (AKR) 

EX ( 1 ) - (0,1.0) * (AKR - AKL) * Z 

EX (2) - (0, -1.0) * (AKR - AKL) * Z 

EX (3) - (0, -1.0) * (AKR + AKL) * Z 

EX (4) - (0,1.0) * (AKR + AKL) * Z 

RAT - ML / MR 
RAT - SQRT(RAT) 

T2(l) - 0.5 * RAT * (1 + AKR/ AKL) * EXP (EX (1)) 

T2 (2) - 0.5 * RAT * (1 + AKR/ AKL) * EXP (EX (2)) 

T2 (3) - 0.5 * RAT * (1 - AKR/ AKL) * EXP (EX (3)) 

T2 ( 4 ) - 0.5 * RAT * (1 - AKR/ AKL) * EXP (EX (4)) 

TR ( 1 ) - Tl (1) *T2 (1) + Tl (3) *T2 (4 ) 

TR (2 ) - Tl (4) *T2 (3) + T1(2)*T2(2) 

TR ( 3 ) - Tl (1) *T2 (3) + Tl (3) *T2 (2) 

TR (4 ) = Tl (4) *T2 (1) + Tl(2)*T2(4) 

Tl (1) - TR ( 1 ) 

Tl (2) - TR(2) 

Tl (3) - TR (3) 

Tl (4) - TR (4) 

RETURN 

END 
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PROGRAM ZIGGURAT 

C This program finds the transmission coefficient T for an 

C incident plane wave of energy E on an asymmetric step potential 

C where the potential gradient (V1-V2) is broken into a series of 

C any desired number of steps NST. The initial and final 

C potential values VO and VF are arbitrary. 

C 

C 

IMPLICIT REAL* 8 (A-H,0-Z) 

DIMENSION Tl(4) 

COMPLEX* 16 T1 
FK » 26.2451 

PRINT*, 1 ENTER VI IN eV (VI is the barrier height) 1 
READ*, VI 

PRINT*, 'ENTER TOTAL BARRIER WIDTH BW IN nm* 

READ*, BW 

PRINT*, 'ENTER THE APPLIED FIELD STRENGTH F IN eV/nm' 

READ*, F 

PRINT*, 'ENTER DESIRED NUMBER OF STEPS NST' 

READ*, NST 
ZO - 0 
VO - 0 

DLTAV = F * BW 
V2 ■ VI - DLTAV 
VF = VO - DLTAV 
DLTAV * DLTAV / NST 
DLTAZ = BW / NST 
TYPE*, DLTAZ 
PRINT 10, 'E(eV) 1 , f T 1 

10 FORMAT ( / / , IX, A7, IX, A5, IX, / IX, 20( , = f ) ) 

OPEN (UNIT - 16, FILE = ' ZIGGURAT . DAT 1 , STATUS - 'NEW') 

Here starts the energy loop 

DO 20 NE = 1, 2000 
EDIV = VI / 2000 
Z * Z0 


E = 1.5 

* 

VI - 

(NE * EDIV) 

T1 ( 1 ) = 

< 

1.0, 

0.0 ) 

T1 (2 ) = 

< 

1.0, 

0.0 ) 

Tl(3) = 

( 

0.0, 

0.0 ) 

T1 (4 ) - 

< 

0.0, 

0.0 ) 

CALL STPTRN (E, 

VO, VI, FK, Z,T1) 


c This takes care of the leading edge of the barrier. Now a loop 

c inside the barrier will handle the potential drop V1-V2. 

c 

IF (DLTAV .EQ. 0) THEN 
Z - ZO + BW 
VL = V2 
VR = VF 

CALL STPTRN (E, VL, VR, FK, Z, T1 ) 

GO TO 29 
END IF 

DO 1*1, (NST - 1) 

Z = Z + DLTAZ 

VL - VI - (I - 1) * DLTAV 

VR = VL - DLTAV 
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CALL STPTRN(E, VL,VR,FK, Z,T1) 

END DO 
C 

c Now the transfer matrix through the trailing edge of the barrier 

c is found and multiplied by the accumulated T-matrix. 

c 

28 VL = VR 
VR = VF 

CALL STPTRN(E,VL,VR,FK, Z, Tl) 

C 

C the control passes to statement 29 from the block IF 

C above which handles the case of zero-field 

C 

29 T = ABS (Tl (1) ) **2 
T = 1 . 0 / T 
PRINT 30, E,T 

30 FORMAT (F7.2,E12.3) 

WRITE (16,300) E, T 

300 FORMAT (2E15.7) 

20 CONTINUE 

21 END 
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PROGRAM ZIGGEFFMASS 
C 

C This program finds the transmission coefficient T for an 

C incident plane wave of energy E on an asymmetric step potential 

C where the potential gradient (V1-V2) is broken into a series of 

C any desired number of steps NST. The applied field F is 

C arbitrary and determines the barrier shape. The effective mass 

C takes an abrupt jump at the edges of the barrier from mx to my, 

C determined by the aluminum contents x and y put in by the 

C operator. This program has the output file ZIGGEFFMASS.DAT. 

C 

C 

IMPLICIT REAL* 8 (A-H,0-Z) 

REAL* 8 Ml, M2 
DIMENSION T1 (4) 

COMPLEX* 16 T1 
FK = 26.2451 

PRINT*, 'ENTER BARRIER HEIGHT VI IN eV' 

READ*, VI 

PRINT*, 'ENTER THE APPLIED FIELD STRENGTH F IN eV/nm' 

READ*, F 

PRINT*, 'ENTER TOTAL BARRIER WIDTH BW IN nm' 

READ*, BW 

PRINT*, 'ENTER DESIRED NUMBER OF STEPS NST' 

READ*, NST 

PRINT*, 'ENTER A1 CONTENT X OF NARROW GAP MATERIAL' 

READ*, X 

PRINT*, 'ENTER Al CONTENT Y OF WIDE GAP MATERIAL' 

READ*, Y 
ZO = 0 
VO = 0 

VAP = F * BW 
V2 = VI - VAP 

VF = VO - VAP 

Ml = 0.0636 + 0.0552 * X + 0.0092 * X**2 

M2 = 0.0636 + 0.0552 * Y + 0.0092 * Y**2 

PRINT*, 'THIS FIELD RESULTS IN V2=',V2,'AND VF= ' , VF 
DLTAV = (VI - V2 ) 

IF (NST.EQ.l) GO TO 6 
DLTAV = DLTAV/ (NST-1) 

6 TYPE*, DLTAV 

DLTAZ = BW/NST 
TYPE*, DLTAZ 

OPEN (UNIT = 18, FILE = 'ZIGGEFFMASS.DAT', STATUS - 'NEW') 

PRINT 10, 'E(eV) ' , 'T' 

10 FORMAT ( //,1X, A7, IX, A5, IX, / IX, 20('=') ) 

Here starts the energy loop 

DO 20 NE = 1, 2000 
EDIV = Vl/2000 
Z = Z0 


E = 1.5 

* 

VI - 

(NE 

* EDIV) 

Tl ( 1 ) - 

< 

1*0, 

0.0 

) 

Tl (2 ) - 

( 

1.0, 

0.0 

) 

Tl (3) = 

( 

0.0, 

0.0 

) 

Tl (4) - 

( 

0.0, 

0.0 

) 


CALL EFFMASSTEP (E, V0, Vl, Ml, M2, FK, Z, Tl) 
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Z - Z + DLTAZ 
c 

c This takes care of the leading edge of the barrier. Now a loop 

c inside the barrier will handle the potential drop VI -V2 . 

c 

IF (DLTAV .EQ. 0) THEN 
Z - Z0 + BW 
VL * V2 
VR = VF 

CALL EFFMASSTEP (E, VL, VR, M2 , Ml , FK, Z, T1 ) 

GO TO 29 
END IF 

DO 25 VI - VI, V2+DLTAV,- (DLTAV) 

VL - VI 

VR - VI - DLTAV 

CALL EFFMASSTEP (E, VI , VR, M2 , M2 , FK, Z,T1) 

Z * Z + DLTAZ 
25 CONTINUE 

c 

c Now the transfer matrix through the trailing edge of the barrier 

c is found and multiplied by the accumulated T-matrix. 

c 

2 8 VL - V2 

VR = VF 

CALL EFFMASSTEP (E, VL, VR, M2 , Ml , FK, Z , T1 ) 

C 

C the control passes to statement 29 from the block IF 

C above which handles the case of zero-field 

C 

29 T - ABS <T1 (1) ) **2 
T - 1.0/T 

PRINT 30, E,T 

30 FORMAT (F7 . 2, E12 . 3) 

WRITE (18,300) E, T 

300 FORMAT (2E15.6) 

20 CONTINUE 

21 END 


138 



PROGRAM MSB 

C The program ALMSB finds the transmission coefficient T as 

C a function of incident plane wave energy E for multiple 

C step potential barriers of height VI, width BW, separation 

C DSEP , and total number NTEETH for AlGaAs-GaAs. Applied field 

C F = 0. Effective mass varies in wells and barriers. Data are 

C written into output file MSB. DAT 

C 

IMPLICIT REAL* 8 (A-H,0-Z) 

DIMENSION T1 (4 ) 

COMPLEX* 16 T1 
FK = 26.2451 
VO = 0.0 
Z0 = 0.0 

PRINT*, 'ENTER BARRIER HEIGHT VI IN eV' 

READ*, VI 

PRINT*, 'ENTER BARRIER WIDTH, BW, AND SEPARATION DSEP IN nm. ' 
READ*, BW, DSEP 

PRINT*, 'HOW MANY BARRIERS ARE THERE ? (NTEETH) ' 

READ*, NTEETH 
PRINT 10, 'E (eV) * , 'T' 

10 FORMAT (//, IX, A7, 1X,A5, IX, /, IX, 20 ( *=' ) ) 

OPEN ( UNIT = 11, FILE = 'ALMSB.DAT' , STATUS ='NEW') 

DO 20 NE = 1,10000 
EDIV = VI / 10000 
Z = Z0 

E = (1.10* VI) - (NE* EDIV) 

T1 ( 1 ) = (1.0, 0.0) 

T1 (2 ) = (1.0, 0.0) 

T1 (3) - (0.0, 0.0) 

T1 (4) = (0.0, 0.0) 

PERIOD = BW + DSEP 
DO 30 NST - 1, NTEETH 

CALL STPTRN (E,V0,V1, FK, Z, T1 ) 

Z = Z + BW 

CALL STPTRN (E, VI, V0, FK,Z,T1) 

Z = NST * PERIOD + Z0 
30 CONTINUE 

C 

C END TEETH LOOP 

C 

T » ABS (T1 (1) ) **2 
T = 1.0/T 
T = LOGIO(T) 

WRITE (11,300) E, T 


300 

FORMAT (2E15.7) 
PRINT 35, E, T 

35 

FORMAT ( 2E12.3 

20 

CONTINUE 

C 


C 

END OF E-LOOP 

c 

END 
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PROGRAM ALMSB 

C The program ALMSB finds the transmission coefficient T as 

C a function of incident plane wave energy E for multiple 

C step potential barriers of height VI, width BW, separation 

C DSEP, and total number NTEETH for AlGaAs-GaAs. Applied field 

C F - 0. Effective mass varies in wells and barriers. 

C 

IMPLICIT REAL* 8 (A-H,0-Z) 

REAL* 8 MX, MY 
DIMENSION T1 (4) 

COMPLEX* 16 T1 
FK = 26.2451 
VO = 0.0 
Z0 = 0.0 

PRINT*, 'ENTER BARRIER HEIGHT VI IN eV' 

READ*, VI 

PRINT*, 'ENTER BARRIER WIDTH, BW, AND SEPARATION DSEP IN nm. * 
READ*, BW, DSEP 

PRINT*, 'HOW MANY BARRIERS ARE THERE ? (NTEETH)' 

READ*, NTEETH 

PRINT*, 'Al CONCENTRATION IN BARRIERS, X,?' 

READ*, X 

PRINT*, 'Al CONCENTRATION IN WELLS, Y,?' 

READ*, Y 

MX = 0.0636 + 0.0552 * X + (0.0092 * X**2) 

MY = 0.0636 + 0.0552 * Y + (0.0092 * Y**2) 

PRINT 10, 'E (eV) ' , 'T' 

10 FORMAT (//, IX, A7, IX, A5, IX, /, IX, 20 ('=' ) ) 

OPEN ( UNIT - 11, FILE = ' ALMSB . DAT ' , STATUS -'NEW') 

DO 20 NE - 1,10000 
EDIV - Vl/10000 
Z - Z0 

E - (1.10 * VI) - ( NE * EDIV) 

T1 ( 1 ) = (1.0, 0.0) 

T1 (2 ) = (1.0, 0.0) 

Tl(3) = (0.0, 0.0) 

Tl (4 ) = (0.0, 0.0) 

PERIOD = BW + DSEP 
DO 30 NST - 1, NTEETH 

CALL EFFMASSTEP (E, V0 , VI , MY, MX, FK, Z, Tl) 

Z = Z + BW 

CALL EFFMASSTEP (E, VI , V0 , MX, MY, FK, Z,T1) 

Z - NST * PERIOD + Z0 


30 

C 

CONTINUE 

C 

END TEETH LOOP 

T - ABS (Tl (1) ) **2 
T = 1.0/T 
T = LOG10 (T) 

WRITE (11,300) E, T 

300 

FORMAT (2E15.7) 
PRINT 35, E, T 

35 

FORMAT ( 2E12.3 ) 

20 

C 

CONTINUE 

C 

END OF E-LOOP 
END 
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PROGRAM SAWTOOTH 


c 

c 

c This program calculates the quantum mechanical transmission 

c coefficient through a series of sawtooth potential barriers 

c under an applied electic field resulting in a voltage drop 

c over the length of the entire barrier of VAP . This program 

c uses the subroutine STPTRN . A data file of T versus E is 

c generated in the output file SAWTOOTH.DAT. 

c 
c 

IMPLICIT REAL* 8 (A-H, 0-Z) 

DIMENSION T1 (4) 

COMPLEX* 16 T1 
FK = 26.2451 
ZOO = 0 
V00 = 0 

PRINT*, 'HOW MANY SAWTEETH, AND HOW WIDE ARE THEIR BASES IN nm ?' 
READ*, NTH, BSE 

PRINT*, 'HOW HIGH IS THE TRAILING EDGE OF EACH TOOTH (eV)?' 

READ*, VI 0 
EDIV = V10 / 2000 

PRINT*, 'ENTER THE APPLIED FIELD IN eV/nm' 

READ*, F 

VAP = F * BSE * NTH 

IF (VAP .EQ. VI 0) VAP = (VAP + EDIV / 2) 

NST = 50 + NINT (ABS ( (- (VAP / NTH) + V10) / 0.01)) 

DLTAZ = BSE / NST 

DLTAV = (-(VAP / NTH) + V10) / NST 

PRINT 10, 'E(eV) ' , 'T' 

10 FORMAT (//, IX, A7 , IX, A5, IX, /, IX, 20 ( * = * ) ) 

OPEN (UNIT = 55, FILE = 'SAWTOOTH.DAT', STATUS = 'NEW') 

DO NE = 1, 2000 

E = 1.10 * V10 - (NE * EDIV) 

T1 ( 1 ) = (1.0, 0.0) 

Tl(2) = (1.0, 0.0) 

T1 ( 3 ) - (0.0, 0.0) 

T1 ( 4 ) = (0.0, 0.0) 

DO NT = 1, NTH 

V0 = V00 - (VAP / NTH) * (NT - 1) 

Z = ZOO + (NT - 1) * BSE 
DO I = 1, NST 

VL - V0 + (I - 1) * DLTAV 

VR = VL + DLTAV 

CALL STPTRN (E, VL, VR, FK, Z, Tl) 

Z = Z + DLTAZ 
END DO 
VL - VR 

VR = V00 - (VAP / NTH) * NT 
CALL STPTRN (E, VL, VR, FK, Z, Tl) 

END DO 

T = ABS (Tl (1) ) **2 
T = (1.0 / T) 

T - LOGIO(T) 

PRINT 100, E, T 
100 FORMAT (F7.2, E12.4) 

WRITE (55,200) E,T 
200 FORMAT (2E15.8) 
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END DO 
END 
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PROGRAM ALSTAIR 


c 

c 

c This program calculates the quantum mechanical transmission 

c coefficient through a series of sawtooth potential barriers 

c under an applied electic field resulting in a voltage drop 

c over the length of the entire barrier of VAP. This program 

c uses the subroutine EFFMASSTEP. It takes into 

c account the changes in effective mass of the electron in each 

c layer in the superlattice structure, it is intended for use 

c in the alas-gaas system, (compare to esaki-tsu) 

c 
c 

IMPLICIT REAL* 8 (A-H, 0-Z) 

DIMENSION T1 (4) 

COMPLEX* 16 T1 

REAL* 8 MMIN f MMAX,ML,MR 

FK - 26.2451 

V00 = 0 

ZOO = 0 

PRINT*, 'HOW MANY SAWTEETH, AND HOW WIDE ARE THEIR BASES IN nm ?' 
READ*, NTH, BSE 

PRINT*, 'HOW HIGH IS THE TRAILING EDGE OF EACH TOOTH (eV) ? ' 

READ*, VI 0 
EDIV = V10 / 2000 

PRINT*, 'WHAT IS THE APPLIED FIELD F? (IN eV/nm) ' 

READ*, F 

VAP = F * BSE * NTH 

IF (VAP .EQ. VI 0) VAP = (VAP + EDIV / 2) 

NST = 50 + NINT (ABS ( {- (VAP / NTH) + V10)/0.01)) 

DLTAZ = BSE / NST 

DLTAV = (-(VAP / NTH) + V10) / NST 

PRINT*, 'WHAT IS THE MAX A1 CONCENTRATION' 

PRINT*, 'ANYWHERE IN THE STRUCTURE?' 

READ*, XMAX 

MMAX = 0.0636 + 0.0552 * XMAX + <0.0092 * XMAX**2) 

PRINT*, 'WHAT IS THE MIN A1 CONCENTRATION' 

PRINT*, 'ANYWHERE IN THE STRUCTURE?' 

READ*, XMIN 

MMIN = 0.0636 + 0.0552 * XMIN + (0.0092 * XMIN**2) 

DLTAAL = (XMAX - XMIN) / NST 
PRINT 10, 'E(eV)', 'LOGIO(T)' 

10 FORMAT (//, IX, A7 , IX, A15, IX, /, IX, 24 ('=')) 

OPEN (UNIT - 59, FILE = 'ALSTAIR.DAT' , STATUS = 'NEW') 

DO NE = 1, 2000 

E = 1.1 * VI 0 - NE * EDIV 
T1 ( 1 ) - (1.0, 0.0) 

T1 (2) = (1.0, 0.0) 

T1 (3) - (0.0, 0.0) 

Tl<4) = (0.0, 0.0) 

DO NT = 1, NTH 

V0 = V00 - (VAP / NTH) * (NT - 1) 

Z = ZOO + (NT - 1) * BSE 
DO I = 1, NST 

VL = V0 + (I - 1) * DLTAV 

VR = VL + DLTAV 

XL = XMIN + (I - 1) * DLTAAL 

ML = 0.0636 + 0.0552 * XL + (0.0092 * XL**2) 
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XR = XL + DLTAAL 

MR - 0.0636 + 0.0552 * XR + (0.0092 * XR**2) 
CALL EFFMASSTEP ( E, VL, VR, ML, MR, FK, Z, Tl) 

Z - Z + DLTAZ 
END DO 

Z = ZOO + NT * BSE 
VL = VR 

VR = V00 - (VAP / NTH) * NT 

CALL EFFMASSTEP (E, VL, VR, MMAX, MMIN, FK, Z, Tl) 
END DO 

T = ABS (Tl (1) ) **2 
T = (1.0 / T) 

T = LOG10 (T) 

PRINT 100, E, T 
100 FORMAT (F7.3, E12.4) 

WRITE (59,200) E, T 
200 FORMAT (2E15.8) 

END DO 
END 
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PROGRAM JVMSB 
C 
C 

C This program calculates the quantum mechanical tunneling current 

C through a series of step potential barriers under an applied 

C electric field resulting in a voltage drop of VAP over the length 

C of the array. It assumes uniform effective mass of mo throughout 

C the structure. Data files of LOG10 (Tunneling current) v.s. applied 

C voltage are written into the output file is JVMSB.DAT 

C 
C 

IMPLICIT REAL* 8 (A-H, 0-Z) 

DIMENSION T1 (4) 

COMPLEX* 16 Tl, AKI , AKT 
FK » 26.2451 
V00 = 0.0 
Z0 - 0.0 
E = 0.005 

PRINT*, 1 ENTER THE NUMBER OF BARRIERS 1 
READ*, NBAR 

PRINT*, ’ENTER THE BARRIER HEIGHTS (eV) , WIDTHS (nm) , 

+ AND SEPARATIONS (nm) 1 
READ*, VI 0, BW, DSEP 
NPER = (NBAR - 1) 

PER * BW + DSEP 

PRINT 10, ’ VAP (eV) ' , 1 LOG10 (TUNCUR) 1 

10 FORMAT ( //, A9, IX, A15, IX, / 1X,24(’~*)) 

OPEN (UNIT - 89, FILE = ’JVMSB.DAT 1 , STATUS * ’NEW’) 

C 

C Here’s the voltage loop 
C 

DO NV - 1, 200 
VDIV = 0.01 
VAP = 2.0 - NV * VDIV 
F = VAP / (NPER * PER + BW) 

NST = 25 + IDNINT (F * BW / 0.001) 

NWST - NST * (DSEP / BW) 

DLTAV = (F * BW) / NST 

DLTAZ - BW / NST 

WDLTAV - (F * DSEP) / NWST 

WDLTAZ * DSEP / NWST 

Tl ( 1 ) = (1. 0,0.0) 

Tl (2 ) - (1. 0,0.0) 

Tl (3) - (0.0, 0.0) 

Tl (4 ) - (0.0, 0.0) 

C 

C Here’s the period loop 
C 

DO N = 1, NPER 

Z = Z0 + (N - 1) * PER 

VL = V00 - F * (N - 1) * PER 
VR = V10 - F * (N - 1) * PER 
CALL STPTRN (E, VL, VR, FK, Z, Tl) 

C 

C Here’s the top edge of the barrier loop 
C 

DO I - 1 , (NST - 1) 
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Z - Z + DLTAZ 

VL - V10 - (F * (N - 1) * PER) - (I - 1) * DLTAV 

VR - VL - DLTAV 

CALL STPTRN {E, VL, VR, FK, Z, Tl) 

END DO 

Here's the trailing edge of the barrier 

Z - ZO + ((N - 1) * PER) + BW 

VL = V10 - F * ( (N - 1) * PER + BW) + DLTAV 

VR - V00 - F * ( (N - 1) * PER + BW) 

CALL STPTRN (E, VL, VR, FK, Z, Tl) 

Here's the well loop 

DO J = 1, (NWST - 1) 

Z - Z + WDLTAZ 

VL “ V00 - F * ( (N - 1) * PER + BW) - (J - 1) * WDLTAV 

VR = VL - WDLTAV 

CALL STPTRN (E, VL, VR, FK, Z, Tl) 

END DO 
END DO 

The following are the computations for the last barrier 

Z - NPER * PER 
VL = VR 

VR = V10 - (F * NPER * PER) 

CALL STPTRN (E, VL, VR, FK, Z, Tl) 

DO M = 1, (NST - 1) 

Z - Z + DLTAZ 

VL - V10 - (F * NPER * PER) - (M - 1) * DLTAV 

VR = VL - DLTAV 

CALL STPTRN (E, VL, VR, FK, Z, Tl) 

END DO 

Z - NPER * PER + BW 

VL = V10 - F * {(NPER * PER) + BW) + DLTAV 
VR = V00 - F * ( (NPER * PER) + BW) 

CALL STPTRN (E, VL, VR, FK, Z, Tl) 

This takes care of all the barriers 

AKI = FK * E 
AKI - SQRT(AKI) 

AKT - FK * (E - VR) 

AKT - SQRT(AKT) 

T - ABS (Tl (1) ) **2 
T - 1.0 / T 

TUNCUR = (AKT / AKI) * T 
TUNCUR = LOG10 (TUNCUR) 

PRINT 20, VAP, TUNCUR 
20 FORMAT ( F7.3, E12.4) 

WRITE (89,30) VAP, TUNCUR 
30 FORMAT (2E15.8) 

END DO 
END 
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PROGRAM JVALMSB 


C 

C 

C This program calculates the quantum mechanical tunneling current 

C through a series of step potential barriers under an applied 

C electric field resulting in a voltage drop of VAP over the 

C length of the array. 

C This program uses the subroutine EFFMASSTEP. It is intended for 

C use with GaAs-Al (x) Ga (1-x) As superlattices, taking account of the 
C effective mass of the electron in each layer of the superlattice. 
C Data files of LOGIO (Tunneling current) vs. applied voltage are 
C written into the output file JVALMSB . DAT 

C 
C 

IMPLICIT REAL* 8 (A-H, 0-Z) 

REAL* 8 MX, MY, MI, MT 
DIMENSION T1 (4 ) 

COMPLEX* 16 Tl, AKI, ART 
FK - 26.2451 
V00 = 0.0 
Z0 - 0.0 
E = 0.005 

PRINT*, 'ENTER THE NUMBER OF BARRIERS' 

READ*, NBAR 

PRINT*, 'ENTER THE BARRIER HEIGHTS (eV) , WIDTHS (nm) 

+ AND SEPARATIONS (nm) ' 

READ*, VI 0, BW, DSEP 

PRINT*, 'Remember, the A1 concentration in the barrier is 
+ higher than in the well!!!!' 

PRINT*, 'WHAT IS THE A1 CONCENTRATION IN THE BARRIER, XMAX? ' 

READ*, XMAX 

PRINT*, 'WHAT IS THE A1 CONCENTRATION IN THE WELL, XMIN? ' 

READ*, XMIN 

MX = 0.0636 + 0.0552 * X + 0.0092 * X**2 
MY - 0.0636 + 0.0552 * Y + 0.0092 * Y**2 
MI - MY 
MT = MY 

NPER - (NBAR - 1) 

PER = BW + DSEP 

PRINT 10, 'VAP(eV)', 'LOG10 (TUNCUR) ' 

10 FORMAT ( //,A9, 1X,A15, IX, / 1X,24('~')) 

OPEN (UNIT - 86, FILE = ' JVALMSB . DAT ' , STATUS = 'NEW') 

Here's the voltage loop 

DO NV = 1, 200 
VDIV = 0.01 
VAP = 2.0 - NV * VDIV 
F = VAP / (NPER * PER + BW) 

NST - 25 + IDNINT (F * BW / 0.001) 

NWST = NST * (DSEP / BW) 

DLTAV = (F * BW) / NST 
DLTAZ = BW / NST 
WDLTAV = (F * DSEP) / NWST 
WDLTAZ = DSEP / NWST 
Tl ( 1 ) = (1.0, 0.0) 

Tl (2) = (1.0, 0.0) 
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Tl<3) = (0.0, 0.0) 

Tl (4) = (0.0, 0.0) 

C 

C Here's the period loop 
C 

DO N = 1, NPER 

Z = Z0 + (N - 1) * PER 

VL = V00 - F * (N - 1) * PER 

VR = V10 - F * (N - 1) * PER 

CALL EFFMASSTEP (E, VL, VR, MY, MX, FK, Z, Tl) 

C 

C Here's the top edge of the barrier loop 
C 

DO I - 1, (NST - 1) 

Z = Z + DLTAZ 

VL = V10 - (F * (N - 1) * PER) - (I - 1) * DLTAV 

VR - VL - DLTAV 

CALL EFFMASSTEP (E, VL, VR, MX, MX, FK, Z, Tl) 

END DO 
C 

C Here's the trailing edge of the barrier 
C 

Z = Z0 + ( (N - 1) * PER) + BW 

VL = VI 0 - F * ( (N - 1) * PER + BW) + DLTAV 
VR = V00 - F * ( (N - 1) * PER + BW) 

CALL EFFMASSTEP (E, VL, VR, MX, MY, FK, Z, Tl) 

C 

C Here's the well loop 
C 

DO J = 1, (NWST - 1) 

Z = Z + WDLTAZ 

VL = V00 - F * ( (N - 1 ) * PER + BW) - (J - 1) * WDLTAV 

VR = VL - WDLTAV 

CALL EFFMASSTEP (E, VL, VR, MY, MY, FK, Z, Tl) 

END DO 
END DO 
C 

C The following are the computations for the last barrier 
C 

Z = NPER * PER 
VL - VR 

VR = VI 0 - (F * NPER * PER) 

CALL EFFMASSTEP (E, VL, VR, MY, MX, FK, Z, Tl) 

DO M = 1, (NST - 1) 

Z = Z + DLTAZ 

VL = VI 0 - (F * NPER * PER) - (M - 1) * DLTAV 

VR = VL - DLTAV 

CALL EFFMASSTEP (E, VL, VR, MX, MX, FK, Z, Tl) 

END DO 

Z = NPER * PER + BW 

VL = V10 - F * ((NPER * PER) + BW) + DLTAV 
VR = V00 - F * ( (NPER * PER) + BW) 

CALL EFFMASSTEP (E, VL, VR, MX, MY, FK, Z, Tl) 

C 

C This takes care of all the barriers 
C 

AKI = FK * MI * E 
AKI = SQRT(AKI) 
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AKT = FK * MT * (E - VR) 

ART - SQRT(AKT) 

T - ABS (T1 (1) ) **2 
T - 1 . 0 / T 

TUNCUR = (AKT / AKI) * (MI / MT) * T 
TUNCUR = LOGIO (TUNCUR) 

PRINT 20, VAP, TUNCUR 
20 FORMAT ( F7.3, E12.4) 

WRITE (86,30) VAP, TUNCUR 
30 FORMAT (2E15.8) 

END DO 
END 
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PROGRAM JVSAWTOOTH 


C 

c 

c This program calculates the quantum mechanical tunneling 

c current through a series of sawtooth potential barriers 

c under an applied voltage VAP. It uses the subroutine STPTRN. 

c Data files of LOG10 (tunneling current) versus applied voltage 

c VAP are created under the filename JVSAWTOOTH.DAT 

c 
c 

IMPLICIT REAL* 8 (A-H, 0-Z) 

DIMENSION T1 (4) 

COMPLEX* 16 Tl, AKI , AKT 
FK - 26.2451 
V00 = 0 
ZOO - 0 
E * 0.005 

PRINT* , ’ HOW MANY SAWTEETH, AND HOW WIDE ARE THEIR BASES IN nm ? f 

READ* , NTH, BSE 

PRINT*, ? HOW HIGH IS THE TRAILING EDGE OF EACH TOOTH (eV) ? 1 
READ*, VI 0 

THIS TAKES CARE OF THE INPUT STATEMENTS 

PRINT 10, ’VAP (V)’, ’LOGIO (TUNCUR) ’ 

10 FORMAT (//, IX, A9, IX, A15, IX, /, IX, 24 (’-’)) 

OPEN (UNIT - 51, FILE * 'JVSAWTOOTH.DAT’, STATUS - ’NEW’) 

DO NV « 1, 200 
VDIV = 0.01 
VAP * 2.0 - VDIV * NV* 

IF (VAP .EQ. VI 0) VAP = VAP + VDIV / 2 

NST « 25 + NINT (ABS ( (- (VAP / NTH) + V10) / 0.001))' 

DLTAZ * BSE / NST 

DLTAV - ( - (VAP / NTH) + VI 0) / NST 

Tl ( 1 ) * (1.0, 0.0) 

Tl (2) * (1.0, 0.0) 

Tl (3) = (0.0, 0.0) 

Tl (4 ) - (0.0, 0.0) 

DO NT - 1, NTH 

V0 - V00 - (VAP / NTH) * (NT - 1) 

Z « ZOO + (NT - 1) * BSE 
DO 1*1, NST 

VL - V0 + (I - 1) * DLTAV 
VR - VL + DLTAV 

CALL STPTRN (E, VL, VR, FK, Z, Tl) 

Z ~ Z + DLTAZ 
END DO 

Z « ZOO + NT * BSE 
VL = VR 

VR « V00 - (VAP / NTH) * NT 
CALL STPTRN (E, VL, VR, FK, Z, Tl) 

END DO 

AKI * FK * E 
AKI - SQRT(AKI) 

AKT - FK * (E - VR) 

AKT = SQRT(AKT) 

T - ABS (Tl (1) ) **2 
T - (1.0 / T) 
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TUNCUR - (AKT / AKI ) * T 
TUNCUR = LOG 10 (TUNCUR) 
PRINT 100, VAP, TUNCUR 
100 FORMAT (F7.3, E12.4) 

WRITE (51,200) VAP, TUNCUR 
200 FORMAT (2E15.8) 

END DO 
END 
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PROGRAM JVALSTAIR 


c This program calculates the quantum mechanical tunneling 

c current through a series of sawtooth potential barriers 

c under an applied electic field resulting in a voltage drop 

c over the length of the entire barrier of vap. This program 

c uses the subroutine EFFMASSTEP. It takes into 

c account the changes in effective mass of the electron in each 

c layer in the superlattice structure. It is intended for use 

c on GaAs-Al (x) Ga (1-x) As sawtooth superlattices, 

c Data files of LOGIO (tunneling current) versus applied voltage 

c VAP are in output file JVALSTAIR.DAT 

c 

IMPLICIT REAL *8 (A-H, 0-Z) 

REAL* 8 Ml, M2 , ML, MR, MI, MT 
DIMENSION T1 (4 ) 

COMPLEX* 16 Tl, AKI, AKT 
FK = 26.2451 
VOO = 0 
ZOO = 0 
E = 0.005 

PRINT*, 'HOW MANY SAWTEETH, AND HOW WIDE ARE THEIR BASES IN nm ?' 
READ*, NTH, BSE 

PRINT*, 'HOW HIGH IS THE TRAILING EDGE OF EACH TOOTH (eV)?' 

READ*, VI 0 

PRINT*, 'WHAT IS THE MAX A1 CONCENTRATION' 

PRINT*, 'ANYWHERE IN THE STRUCTURE?* 

READ*, XMAX 

PRINT*, 'WHAT IS THE MIN Al CONCENTRATION' 

PRINT*, 'ANYWHERE IN THE STRUCTURE?* 

READ*, XMIN 

Ml = 0.0636 + (0.0552 * XMAX) + (0.0092 * XMAX**2) 

M2 = 0.0636 + (0.0552 * XMIN) + (0.0092 * XMIN**2) 

MI = M2 
MT = M2 

PRINT 10, 'VAP(eV)', 'LOG10 (TUNCUR) ' 

10 FORMAT (//, IX, A9, IX, A15, IX, /, IX, 24 {'—')) 

OPEN (UNIT = 60, FILE = 'JVALSTAIR.DAT*, STATUS = 'NEW') 

DO NV = 1, 200 
VDIV =0.01 
VAP = 2.0 - VDIV * NV 

IF (VAP .EQ. VI 0) VAP = (VAP + VDIV / 2) 

NST = 25 + NINT (ABS ( (- (VAP / NTH) + V10)/0.001)) 

DLTAZ = BSE / NST 

DLTAV = (-(VAP / NTH) + V10) / NST 

DLTAAL = (XMAX - XMIN) / NST 
Tl ( 1 ) = (1.0, 0.0) 

Tl (2 ) = (1.0, 0.0) 

Tl (3) = (0.0, 0.0) 

Tl (4 ) = (0.0, 0.0) 

DO NT = 1, NTH 

V0 = VOO - (VAP / NTH) * (NT - 1) 

Z = ZOO + (NT - 1) * BSE 


DO I 

= 1, NST 


VL 

= V0 + (I - 1) * 

DLTAV 

VR 

= VL + DLTAV 


XL 

= XMIN + (I - 1) 

* DLTAAL 
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ML = 0.0636 + 0.0552 * XL + (0.0092 * XL**2) 
XR - XL + DLTAAL 

MR - 0.0636 + 0.0552 * XR + (0.0092 * XR**2) 
CALL EFFMASSTEP ( E, VL, VR, ML, MR, FK, Z, Tl) 
Z = Z + DLTAZ 
END DO 

Z = ZOO + NT * BSE 
VL - VR 

VR - V00 - (VAP / NTH) * NT 

CALL EFFMASSTEP (E, VL, VR, MR, M2, FK, Z, Tl) 
END DO 

AKI = FK * MI * E 
AKI - SQRT(AKI) 

AKT = FK * MT * (E - VR) 

AKT -= SQRT(AKT) 

T - ABS(T1(1) )**2 
T - (1.0 / T) 

TUNCUR - (AKT / AKI) * (MI / MT) * T 
TUNCUR - LOG 10 (TUNCUR) 

PRINT 100, VAP, TUNCUR 
FORMAT (F7.3, E12.4) 

WRITE (60,200) VAP, TUNCUR 
FORMAT (2E15.8) 

END DO 
END 
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PROGRAM WELL 


C 

C This program finds the total transfer matrix element Mil 

C for a single quantum well of depth VI and width a, as a 

C function of the energy below the top of the well. A constant 

C effective mass m* = mo is assumed. Data files of LOGlO(Mll) 

C versus E are generated. 

C 

IMPLICIT REAL* 8 (a-h,o-z) 

REAL* 8 Mil 
DIMENSION T1 (4) 

COMPLEX* 16 T1 
FK - 26.2451 
ZO = 0 
VO = 0 

PRINT*, 'ENTER WELL DEPTH VI IN eV (ENTER ABS VALUE) ' 

READ*, VI 

PRINT*, 'ENTER WELL WIDTH, a, IN NANOMETERS' 

READ*, A 

OPEN (UNIT - 80, FILE - 'WELL', STATUS = 'NEW') 

PRINT 10, 'E(eV) ', 'LOGlO(Mll) ' 

10 FORMAT (//, IX, A7, IX, All, IX, /,1X,24 ('=') ) 

VI = -VI 

EDIV = VI / 2000 
DO 20 NE - 1,2000 
E - NE * EDIV 
Z - Z0 

T1 (1) - (1.0, 0.0) 

T1 (2) = (1.0, 0.0) 

T1 (3) - (0.0, 0.0) 

T1 (4 ) - (0.0, 0.0) 

CALL STPTRN (E, V0, VI, FK, Z, Tl) 

Z - Z + A 

CALL STPTRN (E, V1,V0, FK, Z, Tl) 

Mil = ABS (Tl (1) ) 

Mil = LOG10 (Mil ) 

WRITE (80,300) E, Mil 
300 FORMAT (2E15.6) 

PRINT 35, E, Mil 
35 FORMAT (F7.2, E12.3) 

20 CONTINUE 

END OF E-LOOP 

END 
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PROGRAM ALGAASWELL 


"• This program calculates the total transfer matrix element 

: Mil for a single A1 (x) Ga (1-x) As-Al (y) Ga (1-y) As quantum well 

: of depth VI and width a, as a function of the energy below the 

: top of the well. The input well depth must be in accord with 

2 the values used for x and y, where x and y are the aluminum 

Z concentrations outside and inside the well, respectively . Data 

2 files of LOGl 0 (Mil ) versus E are generated. 

IMPLICIT REAL* 8 (A-H,0-Z) 

REAL* 8 Ml, M2, Mil 
DIMENSION T1 (4) 

COMPLEX* 16 T1 
FK = 26.2451 
ZO = 0 
VO = 0 

PRINT*, 'ENTER WELL DEPTH VI IN eV (ENTER ABS VALUE)' 

READ*, VI 

PRINT*, 'ENTER WELL WIDTH, a, IN NANOMETERS' 

READ * A 

PRINT*, 'ENTER MIN AL CONTENT, y, (IN THE WELL)' 

READ*, XMIN 

PRINT*, 'ENTER MAX AL CONTENT, x, (OUTSIDE THE WELL)' 

READ*, XMAX 

Ml = 0.0636 + (0.0552 * XMAX) + (0.0092 * XMAX**2) 

M2 = 0.0636 + (0.0552 * XMIN) + (0.0092 * XMIN**2) 

OPEN (UNIT = 73, FILE = 'ALGAASWELL', STATUS - 'NEW') 

PRINT 10, ' E (eV) ' , ' LOG10 (Mil) ' 

10 FORMAT (//,1X,A7,1X,A11,1X,/ ,1X,24 ('=') ) 

VI = -VI 

EDIV = VI / 2000 
DO 20 NE = 1,2000 
E - NE * EDIV 
Z = Z0 

T1 ( 1 ) = (1.0, 0.0) 

T1 (2) = (1.0, 0.0) 

Tl (3) -= (0.0, 0.0) 

T1 (4) = (0.0, 0.0) 

CALL EFFMASSTEP (E, V0 , Vl , Ml , M2, FK, Z,T1) 

Z = Z + A 

CALL EFFMASSTEP (E, VI , V0 , M2 , Ml , FK, Z,T1) 

Mil = ABS (Tl (1) ) 

Mil = LOG10 (Mil) 

WRITE (73,300) E, Mil 
300 FORMAT (2E15.6) 

PRINT 35, E, Mil 
35 FORMAT (F7.2, E12.3) 

20 CONTINUE 

END OF E-LOOP 

END 
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PROGRAM STARKMASS 


C 

C This program calculates the total transfer matrix element 

C Mil in a single A1 (x) Ga (1-x) As-Al (y) Ga (1-y) As quantum well, 

C of depth VI and width a as a function of the energy below the 

C top of the well and the applied field F. This program assumes 

C the applied field is restricted to the area of the well. 

C The input value of VI must be compatible with the values of 

C x and y, the aluminum concentrations outside and inside the 

C well, respectively. 75 intervals are used. A data file of 

C LOGlO(Mll) versus E is generated. 

C 

C 

IMPLICIT REAL* 8 (A-H,0-Z) 

REAL* 8 Ml, M2, Mil 
DIMENSION Tl<4) 

COMPLEX* 16 Tl 
FK - 26.2451 
ZO - 0 
VO - 0 
NST - 75 

PRINT*, 1 ENTER WELL DEPTH | VI | IN eV 1 
READ*, VI 

PRINT*, 'ENTER WELL WIDTH a IN nm' 

READ*, A 

PRINT*, 'ENTER THE APPLIED FIELD STRENGTH F IN eV/nm' 

READ*, F 

PRINT*, 'ENTER A1 CONTENT x OUTSIDE WELL' 

READ*, XMAX 

PRINT*, 'ENTER Al CONTENT y INSIDE WELL' 

READ*, XMIN 

Ml = 0.0636 + 0.0552 * XMAX + (0.0092 * XMAX**2) 

M2 = 0.0636 + 0.0552 * XMIN + (0.0092 * XMIN**2) 

VI - -VI 

VAP - F * A 

V2 = VI - VAP 

VF = V0 - VAP 

DLTAV = ABS (VI - V2 ) 

DLTAV - DLTAV/ (NST) 

DLTAZ - A/NST 

OPEN (UNIT =13, FILE = ' STARKMASS . DAT ' , STATUS = 'NEW') 

PRINT 10, 'E (eV) ' , 'MU' 

10 FORMAT ( //,1X, A7 , IX, All, IX, / IX, 24('- f ) ) 

C 

C Here starts the energy loop 

C 

EDIV = Vl/2000 
DO 20 NE = 1, 2000 
E = (NE * EDIV) 

Z = Z0 


Tl(l) - 

( 

1.0, 

0.0 

) 

Tl (2 ) - 

( 

1.0, 

0.0 

) 

Tl (3) - 

( 

0.0, 

0.0 

) 

Tl (4 ) - 

( 

0.0, 

0.0 

) 


CALL EFFMASSTEP (E, V0 , VI , Ml , M2 , FK, Z, Tl) 
IF (DLTAV .EQ. 0) THEN 
Z = Z0 + A 
VL - V2 


157 



non 


vr = VF 

CALL EFFMASSTEP (E, VL, VR, M2 , Ml , FK, Z, T1 ) 

GO TO 29 
END IF 

DO 1=1, (NST - 1) 

Z - Z + DLTAZ 

VL - VI - (I - 1) * DLTAV 

VR = VL - DLTAV 

CALL EFFMASSTEP (E, Vl , VR, M2 , M2 , FK, Z, T1 ) 

END DO 

28 VL = V2 
VR = VF 

Z - ZO + A 

CALL EFFMASSTEP <E, VL, VR, M2 , Ml , FK, Z, T1 ) 

C 

C the control passes to statement 29 from the block IF 

C above which handles the case of zero-field 

C 

29 Mil = ABS (T1 (1) ) 

Mil - LOGIO (Mil ) 

PRINT 30, E, Ml 1 

30 FORMAT (F7.3,E12.3) 

WRITE (13,300) E, Mil 

300 FORMAT (2E15.7) 

20 CONTINUE 

END OF E-LOOP 

21 END 
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o o o n n 


PROGRAM AIRYWELL 


C 

C THIS PROGRAM CALCULATES THE RATIO OF COEFFICIENTS OF THE WAVE 
C FUNCTION ON THE DOWNSTREAM SIDE OF A SINGLE QUANTUM WELL, (F/E) , 

C AS A FUNCTION OF INCIDENT ENERGY AND APPLIED ELECTRIC FIELD. 

C A PLOT IS THEN GENERATED GIVING THIS RATIO VERSUS ENERGY. WHERE- 

C EVER THE RATIO IS MAXIMIZED, A VIRTUAL STATE EXISTS. THIS PROGRAM 

C ASSUMES A GaAs-Al (0 . 5) Ga (0 . 5) As 3 run WIDE WELL, 0.5 eV DEEP, WITH 
C CONSTANT EFFECTIVE MASS m* = 0.0636. 

C 

C 

DIMENSION EARR (1000) , RARR(1000) 

W = 1.5 
V0 = 0.4 

PRINT*, 'WHAT APPLIED FIELD STRENGTH IN V/nm? ' 

READ*, FAPP 
EMAX = FAPP * W 
EMIN = -V0 - FAPP * W 
FAPP = FAPP * 1.0E9 

OPEN GKS 

CALL GOPKS (6,1) 

OPEN THE WORKSTATION 
CALL GOPWK (1,2,1) 

ACTIVATE THE WORKSTATION 
CALL GACWK (1) 

DO I = 1, 1000 
A « 1 . 2103E-3 
B = 1.2147E6 

EARR ( I ) = EMIN + (1/1000.0) * (EMAX - EMIN) 

X1NEG - A* (-W) * FAPP** (1 . /3. ) + B*EARR ( I ) /FAPP* * (2 . /3 . ) 

X1POS = A* (-W) *FAPP** (1./3. ) + B* (EARR (I) + V0) /FAPP** (2 . /3 . ) 
X2NEG = A* (W) *FAPP** (1 . /3. ) + B* (EARR (I ) + V0 ) /FAPP* * (2 . /3 . ) 
X2POS = A* (W) *FAPP** (1 . /3. ) + B*EARR(I) /FAPP** (2. /3. ) 

C - AI (-X1NEG) *BID (-X1POS) - AID (-X1NEG) *BI (-X1POS) 

D = AI (-X1POS) *AID (-X1NEG) - AID (-X1POS) *AI (-X1NEG) 

P = AI (-X2NEG) /AI (-X2POS) 

Q = BI (-X2NEG) /AI (-X2POS) 

R - BI (-X2POS) /AI (-X2POS) 

S - AID (-X2NEG) /AID (-X2POS) 

T = BID (-X2NEG) /AID (-X2POS) 

U - BID (-X2POS) /AID (-X2POS) 

F - C/D * (S - P) + T - Q 

E » U * (C/D * P + Q) - R * (C/D * S + T) 

RARR(I) ■= F/E 
END DO 

CALL EZXY (EARR, RARR, 1000, ' f/e' ) 

DEACTIVATE AND CLOSE THE WORKSTATION 
CALL GDAWK (1) 

CALL GCLWK (1) 

C CLOSE GKS 

CALL GCLKS 

STOP 

END 
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